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ABSTRACT 

We have developed a new method, K2, optimized for the detection of galaxy clusters in multicolor 
images. Based on the Red Sequence approach, K2 detects clusters using simultaneous enhancements in 
both colors and position. The detection significance is robustly determined through extensive Monte- 
Carlo simulations and through comparison with available cluster catalogs based on two different 
optical methods, and also on X-ray data. K2 also provides quantitative estimates of the candidate 
clusters' richness and photometric redshifts. Initially K2 was applied to the two color (gri) 161 
deg 2 images of the Canada-France-Hawaii Telescope Legacy Survey Wide (CFHTLS- W) data. Our 
simulations show that the false detection rate for these data, at our selected threshold, is only ~ 1%, 
and that the cluster catalogs are ~ 80% complete up to a redshift of z = 0.6 for Fornax-like and richer 
clusters and to z ~ 0.3 for poorer clusters. Based on the g, r and i-band photometric catalogs of the 
Terapix T05 release, 35 clusters/ deg 2 are detected, with 1-2 Fornax- like or richer clusters every two 
square degrees. Catalogs containing data for 6144 galaxy clusters have been prepared, of which 239 
are rich clusters. These clusters, especially the latter, are being searched for gravitational lenses - one 
of our chief motivations for cluster detection in CFHTLS. The K2 method can be easily extended to 
use additional color information and thus improve overall cluster detection to higher redshifts. The 
complete set of K2 cluster catalogs, along with the supplementary catalogs for the member galaxies, 
are available on request from the authors. 

Subject headings: galaxies: clusters: general — methods: miscellaneous — catalogs 



1. INTRODUCTION 

Galaxy clusters are the most massive virialized struc- 
tures in the present day universe, with masses ranging 
from a few times 10 14 M for poor clusters to rich clus- 
ters at mass ~ 10 15 M Q or higher. Occurring at the inter- 
section of mass filaments, clusters continue to build up 
mass thorough accretion of infalling galaxies and inter- 
cluster gas. Since clusters occur at the high end of halo 
mass function, their observed mass function and its evo- 
lution with redshift have therefore found extensive use as 
test grounds for cosmogony and large s cale structure for- 
mation models dVikhlinin et al.l 120091 iMalinovskv et al.l 
120081 iMantz et alj 120081 iPierpaoli et alj 120031) and for 
constraining co smolo gical parameters (IWik et all 120081. 
Gladders etail 120071 Marker et al.ll2007l iSefusatti et all 



2007). In addition, clusters also provide excellent test 



beds for models describing environmental effects on 
galaxy evolution. Since the deep gravitational potential 
well of a cluster retains all the accreted baryonic mat- 
ter, detailed observations of the member galaxy popula- 
tion and of the intra-cluster medium provide important 
clues toward the star formation history and the pres- 
ence and strength of associated feedback processes on 
the evolution of galaxies in high density environments 
(jNagai et all 127)071 I\^itlf2005h . 

Our focus on clusters is for yet another application 
- their high projected surface mass de nsities result in 
large gravitational lens cross sections (jHennawi et alJ 
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l2007LlSmith et alJ l2001). thus making them very efficient 
in lensing high redshift objects which lie in th e direc - 
tion of thei r lines of sight (e.g., lEbeling et alj {2009), 
iSmail et alj (120071) ). Spectac ular cluster lenses, such as 
Abell 2218 (jPello et al.lll992l ). are well known examples 
of this phenomenon. The magnification boost offered 
by these cluster lenses, when combined with modern ob- 
servational techniques offers a viable technique to re- 
solve sub-galactic scales in the high redshift universe; for 
example, Gemini GMOS 2D in t egral field (IFU) spec- 
troscopy dSwinbank et~aT1 l200l l2006h and K eck LRIS 
longslit spectroscopy ( Pettini et al.ll2000l l2002j ) of lensed 
arcs have been used to map the kinematics, intensity of 
star formation, as well as the effect of feedback processes 
on the evolution of high redshift (z ^> 1) galaxies. The 
number of such investigations is limited mainly by the 
lack of a significant sample of bright, lensed galaxies to 
carry out these studies. 

In order to detect new cluster lenses, we have under- 
taken a search for clusters and groups of galaxies and any 
associated lenses in multi-color photometric data with 
large sky coverage. The search is carried out in two 
steps: first detecting galaxy clusters and groups, (which 
have a higher likelihood of lensing background objects 
due to their mass concentrations), and subsequently car- 
rying out a dedicated search for strong lens images in 
these cluster regions. In this publication, planned as the 
first of two papers, we present details of our cluster de- 
tection methodology, results from the Monte Carlo cali- 
brations for contamination and expected completeness of 
the search, and the cluster catalogs generated by the de- 
tector from 161 deg 2 of imaging from the Canada France 
Hawaii Telescope Legacy Survey, Synoptic Wide compo- 
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nent (referred hence as, CFHTLS-W). The lens search 
scheme and the resulting lens catalogs will be presented 
in a succeeding publication. 

It must be added that the public release of data from 
several large recent surveys has consequently encouraged 
the rapid development of automated searches for strong 
lenses, ranging from quas ars lenses in Sloan Digital Sky 
Survey (SDSS) imaging (|Inada et alj 120081 ). to galaxy- 
galaxy scale lenses, e.g., robotic s e arches in HST-ACS 
imaging archives ([Marshall et all f2 009. B olton et all 
[2008f) . in CFHTLS imaging ( jCabanac et all 120071). and 
visual search in Cosmic Evolution Sur vey (COSMOS) 
fields (jFaure et alj 120081 . Uacksonl I2008J) , and to group 
scale lenses, e.g. in CF HTLS imaging { Limousi n et all 
l2009f) and SDSS imaging (jKubo et al.ll2009D . The strong 
lcnsing candidates in these catalogs are fully comple- 
mented by the search we have undertaken for giant arcs 
in cluster scale lenses in CFHTLS- Wide imaging. 

The recently completed CFHTLS project is a collab- 
oration between the CFHT corporation, le Centre Na- 
tional de la Recherche Scientifique de France, le Com- 
missariat a l'Energie Atomique and the National Re- 
search Council of Canada. Organized with a Pi-less open 
structure, this wide field imaging survey comprises of 
three components: Deep Synoptic, Wide Synoptic and 
the Very Wide surveys, each with different scientific ob- 
jectives and observing strategies. The principal goal of 
the Wide Synoptic part, the data of which we use for the 
cluster search, is to map structure formation in the uni- 
verse using cosmological weak lensing; a more complete 
description of the objectives of all three components as 
well as other pertinent information are provided at the 
CFHT website 3 . 

The CFHTLS observations were carried out at the 
Canada-France-Hawaii 3.6m telescope, using the 1 
square degree wide field imager, MegaPrime. The 
CFHTLS-W aimed to map 195 square degrees in the five- 
filter ugriz set 4 and in four different sky patches (referred 
as Wl, W2, W3 and W4) to a depth of i < 24 mag. In ad- 
dition, an image quality constraint of seeing better than 
0".9 in r-filter was imposed so that the stacked images 
could be used to construct shear maps for weak lensing 
analysis, the principal science goal. For the CFHTLS- 
W, the survey strategy was to observe all the fields first 
with only the gri filters, with follow-up observations car- 
ried out in urz. Table [T] is a summary of the completion 
statistics of the four CFHTLS-W fields at the end of the 
survey period in August 2008; the statistics listed in the 
table have been computed using the observational de- 
tails provided at the Terapix 5 website for each field. The 
complete data processing, stacking of images and their 
astrometric and photometric calibrations are carried out 
at Terapix and the complete set of data products, in- 
cluding the photometric catalogs which we use for our 
cluster detection, is archived and distributed through 
the Canadian Astronomy Data Centre 6 (CADC). Our 
principal objective is to implement a cluster detection 

3 http:/ /www. cfht.hawaii.edu/Science/CFHLS/ 

4 We do not use the 'primes' in the notation of these MegaPrime 
u*g'r'i' z' filters, which closely match the Sloan u'g'r'i'z' filters in 
their characteristics 

5 http://terapix.iap.fr/ 

6 http:/ /wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 



method tuned for 2-color imaging data and with a well 
understood selection function, to build a large sample of 
clusters at a median redshift z ~ 0.6. 

To motivate the development of our cluster detec- 
tor, the paper begins with a discussion on various pub- 
lished methods designed for optical data in Section [2l 
We then introduce our cluster detection method in Sec- 
tion [3l highlighting its relevance and specific advan- 
tages for this particular photometric data. The Monte- 
Carlo simulations we have used to test the selection 
bias, completeness and false detection rate of our detec- 
tion algorithm are described in Section |U As a further 
test of our method, we compare our cluster detections 
from the 4 deg 2 of CFHTLS-Deep imaging with those 
in t hree other publishe d catalogs, the XMM-LSS 7 clus- 
ters (jPierre et a l. 2007), the catalogs from the photomet- 
ric re dshift based Adaptive Filter method (jMazure et al.l 
l2007h as wel l as cluster candida tes from the Matched Fil- 
ter catalog (jOlsen et al.1 120071 ); results from these com- 
parisons are presented in Section O We present our 
cluster catalogs from the 161 deg 2 of CFHTLS-W imag- 
ing, with each candidate categorized into high, medium 
and low quality based on the detection likelihood, in 
Section [H We summarize all our results in Section [7] 
and outline improvements in progr ess and future plans . 
We use Five- Year W MAP values (jDunklev et all 120091 
iKomatsu et alj l2009f ) of the cosmological parameters 
throughout. 

2. OPTICAL DETECTION OF CLUSTERS 

Optically selected cluster catalogs spanning ever in- 
creasing redshift intervals continue to be published using 
innovative detection techniques that have kept pace with 
the advances in optical observational technology. In ad- 
dition, cluster catalogs may also be constructed using 
detections at other wavelengths, e.g., with X-rays emit- 
ted by the hot gas trap ped in the deep grav itational po- 
tential of clusters (e.g. iPacaud et al.1 (|2007D ). as well as 
in microwave wa velengths using the thermal Sunyaev- 
Zel'dovich effect (jLaRoque et al.l[2003f ) . to cite just two; 
however, we restrict our present discussion to only op- 
tical methods. K2 draws upon the strengths of several 
earlier methods while addressing their shortcomings, and 
the descriptions and discussions in this section explain 
our motivations for devising this new method and to ac- 
knowledge existing techniques which we have built into 
its design. Readers not needing this background may 
skip directly to Section [3l 

Optical detection methods identify clusters utilizing 
one or more of the characteristic properties of their 
member galaxies to effectively isolate them from their 
field galaxy counterparts. Amongst these, the enhanced 
density of bright galaxies is the most visible marker 
of the presence of a cluster. Used by both Abell 
and Zwicky i n their pioneering work on v isual clus- 
ter detection (|Abelllll958l iZwickv et~ai1 IT968). this den- 
sity en hancement is used by the automated Counts in 
Cells dLidman fc Petersonl I1996D and Vor noi Tessella- 
tion (jKim et al.l I2002L iRamella et al.ll2001f ) methods to 
effectively isolate candidate clusters in single filter imag- 
ing, However, since these methods rely only on posi- 

7 XMM-Newton (X-ray Multi Mirror) Large Scale Structure sur- 
vey 
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tion information, they are prone to contamination due 
to projection effects from foreground and background ob- 
jects, making them a poor choice for surveys, such as the 
CFHTLS-W, for which multi-filter data are available. 

This enhanced galaxy density contrast between clus- 
ter and field environments is further accentuated by the 
'Morphology - Density', (T — E) relation - there is an 
increasing fraction of early type (E+SO) galaxies in the 
higher density cores of clusters while the late type galax- 
ies are preferentially fo und in lower density, field environ- 
ments (|Dressleri ri980). Recent surveys have confirmed 
this morphologic al segregation in clusters even at red- 
shift 1 and above (|Helsdon fe Pon man 2003, Smit h~et alj 
2005). Further, these surveys have also found an evo- 
lution in the (T — S) relation with redshift, with the 
early type galaxy fraction, f(E+so) = 0.7 being lower at 
z = 1 compared to the value of 0.9 observed in the local 
universe. This evolution in the E+S fraction is consis- 
tent w ith the ' Butcher- Oemler' effect (Butcher fc Oemierl 
11984) in which high redshift clusters are observed to have 
an excess of blue galaxies in their cores, indicative of 
ongoing star formation compar ed to their low redshift 
counterparts (jGerke et alJ 120071 ). In combination, these 
two effects manifest themselves as a higher density of 
(E+S0) galaxies in the cores of clusters compared to the 
field, with the colors of these galaxies evolving toward 
the red with advancing cosmic time. 

Based on this model of a concentration of bright 
(E+S0) gala xies in the cores of clusters, the Matched 
Filter (MF) ([Postman et al.l[l996h applies a filter which 
combines the expected (Schechter) luminosity function of 
cluster members as well as their radial (circularly sym- 
metric) distribution to either a pre-defined grid or to in- 
dividual galaxies (above a magnitude threshold) in imag- 
ing data and detects clusters as enhancements in the re- 
sulting density map. In the Matched Filter method, the 
reliance on an 'ideal' cluster model is both a strength 
and a drawback. By fitting a fiducial cluster model 
to the data, the method derives useful cluster param- 
eters, such as the richness and redshift, as part of the 
detection process. The method is ideal for deep, single 
filter data and can be tuned for different cluster mod- 
els. At the same time however, the dependency on a 
cluster model biases the method against clusters which 
do not fit the a ssumed profile, e.g. a Bautz-Morgan 
Type III cluster ([Bautz fc Morgan1ll970l ). without a cen- 
tral galaxy, may be missed by the MF method. This 
drawback becomes greater at higher redshifts, where the 
higher m erger rates produce s ignificant differences among 
clusters ([Cohn fc White 2005). In addition, this method 
too is prone to projection effects since it fits all galaxies, 
including projected field galaxies that lie within the aper- 
ture, to the assumed cluster model. This contamination 
may be significantly reduced by including available color 
information, a technique that is used by the methods de- 
scribed in the following paragraph, and by our detection 
algorithm as well. 

In addition to the segregation by morphology, 
the (E+S0) galaxy population in a cluster, espe- 
cially those within the central core regions, shows a 
strong correlation between color and absolute mag- 
nitude. This Colo r -Magnitude Relation (CMR), 
(Ide Vaucouleursl 119611 iVisvanathan fc Sandagd 119771. 
iBower et al.l 11992. Stanfor d et al.l 11998. Ko dama et al.l 



[19981 IBlakeslee et all [20031 IDe Lucia et all 120071) 1 

is observed as a t ight ridge, the Red Sequence 
(|Gladders fc Y ee 2000) , which extends over several mag- 
nitudes with small scatter (< O.lmag) in a color- 
magnitude diagram (CMD); recent results using near 
infra-red observations indicate that the CMR may have 
already been in place in proto-clusters as early as z ~ 3 
(|Kodama et al.ll2007[) . There is an emerging consensus 
that the observed tight Red Sequence indicates that the 
stellar populations in cluster ellipticals are homogeneous, 
formed very early (z ^> 1) and have been passively evolv- 
ing since then. The observational evidence clearly shows 
that the CMR is universal, being present even in poor 
clusters and is also identical in zero point and slope - 
within photometri c uncertainties - in any two clusters at 
the same redshift ([Gladders fc Yeell2000l iKodama etHI 
l2007h . In addition, brighter, and hence more massive 
galaxies (assuming a simple mass-to-light ratio), are ob- 
served to have redder colors than less massive ones, which 
therefore leads to a small slope of ~ 0.1 mag/mag in this 
color sequence ([Kodama et al.|[l998h . 

Therefore, searching for clusters in color space by using 
the Red Sequence of the (E+S0) cluster members offers a 
complementary avenue to reduce the dependence on an 
ideal cluster model (based only on the magnitude and 
distribution of the member galaxies) and the consequent 
selection bias. Such a color-ba s ed ap proach was first im- 
plemented by I Gladders fc Yed {2000) in the Cluster Red 
Sequence (CRS) method, which identifies clusters as si- 
multaneous overdensities in a combined (color+position) 
space. The overdensity in color space due to the clus- 
tering of member galaxies along the Red Sequence is 
isolated using a series of overlapping color slices on a 
color-magnitude diagram; the slope of each color slice 
is matched to that of a fiducial Red Sequence in a cho- 
sen set of redshift intervals, while the width of each slice 
is set by the photometric error of the available imaging. 
This color+position approach efficiently isolates clusters, 
obtains their (photometric) redshift as part of the detec- 
tion process and also effectively suppresses foreground 
and background contamination independent of the red- 
shift of the candidate cluster. With the need for color 
information however, such methods are restricted to sur- 
veys in which the imaging is available in two and prefer- 
ably more filters. The CRS has been implemented in 
the ongoing Red Sequence Cluster Survey (RCS), a ded- 
icated cluster survey, which uses a matched pair of R c 
and z' filters to isolate the stron g 4000A break in earl y 
type galaxies at redshifts up to 1 ([Gladders fc Yeel l2005) . 

An alter nate color based a pproach, the Cut and En- 
hance CCE. IGoto et~aT1 {HI)) designed for SDSS imag- 
ing, applies a series of color cuts to isolate subsets of 
galaxies with similar colors from three filter (Sloan gri) 
photometric catalogs; this selection based on color clus- 
tering is augmented with the galaxy number density 
computed within cells on a grid. The color-position 
clustering signals are amplified by using a weighting 
scheme to create densit y maps on which 'Sextractor' 
(|Bertin fc Arnoutill996D is used to detect peaks as can- 
didate clusters. 

In the closely related MaxBCG method (|Hansen et al.l 
2005), also developed to work with SDSS imaging, each 
galaxy is tested for the likelihood of being a Bright Clus- 
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ter Galaxy (BCG). The method uses a combination of 
the number of galaxies that lie within a defined color 
slice of the galaxy being tested (thus checking for a Red 
Sequence) and a fit of the galaxy's properties to an em- 
pirical BCG evolutionary track. The likelihood function 
is tested at a series of redshift slices and maxima in the 
likelihood function are used to locate BCGs and the clus- 
ters in which they occur. Based on a BCG evolutionary 
model, this method obtains the richness and redshift of 
the cluster as part of the fit. It must be pointed out how- 
ever that recent results (|De Lucia fc Blaizotl 120071 iBildl 
120071 ) indicate that BCGs in certain rich clusters have 
ongoing star formation with colors up to ~ lmag bluer 
than the red sequence; such clusters would not be identi- 
fied by the MaxBCG. In addition, like the Matched Filter 
method, the MaxBCG too would miss clusters which lack 
a central galaxy, i. e., clusters classified as Bautz- Morgan 
Type III clusters (jBautz fe Morgan! ll970T ). Finally, the 
dependency on an empirical evolutionary track for all 
BCGs adds a further selection bias. 

Knowledge of the redshift of even a subset of the galax- 
ies in the survey area provides a third dimension and a 
very stringent constraint that may be employed for clus- 
ter detection. By incorporating available redshifts, the 
Adaptive Matched Filter (AMF) retains all the powers 
of th e basic Matched Filter algorithm (Postman et al. 
1996) while significantly reducing or removing noise from 
foreground and background contamination; the method 
also provides better estimates of cluster membership and 
richness. However, the dependency on a cluster model 
still leaves it pron e to miss clusters which do not match 
the input model, (|Kim et al.ll2002f ). The C4 algorithm 
(IMiller et al.ll2005fl . implemented on 2600 deg 2 of SDSS 
data, incorporates position, five filter photometry and 
spectroscopic redshifts to identify clusters at over 90% 
completeness and 95% purity for clusters with mass 
> 10 14 Mq and up to redshift 0.12, as shown by extensive 
Monte-Carlo tests. The main challenge in applying this 
method for detecting clusters at higher redshift is the de- 
pendency on spectroscopic redshifts, which at present is 
expensive in terms of observing time. 

Photometric redshifts (photoz), estimated from multi- 
color imaging, offer a viable alternative to utilize red- 
shift information for cluster detecti on without incurring 
the penalty of long exposure times. iMazure et al.l (|2007l ) 
have established a benchmark for this approach using the 
public photometric redshift catalog for CFHTLS-Deep 1 
field, estimat ed from multi-colo r optical and near infra- 
red imaging (II bert et alJ 12006) . For cluster detection, 
IMazure et alJ (|2007l ) use an adaptive kernel to locate 
density enhancements in i-band imaging in overlapping 
photometric redshift slices of width A z = 0.1. In a com- 
parison with the X MM-LSS cluster catalog for this field 
(| Pierre et~a l. 2007), the method is shown to recover 100% 
of the X-ray detections. Given the overlap of CFHTLS- 
Dl with our search area, we com pare the K2 detec tions 
with those of bot h the XMM-LSS (iPierre et aill2007l ) and 
Adaptive Kernel (|Mazure et al.ll2007l ) catalogs in Section 
[51 and therefore provide further pertinent details of both 
these catalogs in the corresponding subsections. 

In an alternative appro ach using photometric redshifts, 
iGavazzi fe Sou cail (2007) combine weak lensing recon- 
structions in the i-band images with the available pho- 



toz catalogs for the four CFHTLS-Deep fields to identify 
clusters as peaks in the convergence maps. Comparison 
with the overlapping XMM-LSS catalogs indicates that 
nine of the fourteen X-ray clusters are secure detections 
with redshifts and velocity disp ersions being measu red as 
part of the detection process. iBerge et~a l. (2008) use a 
similar approach but use wavelets for reconstructing the 
weak lensing maps; these maps are then combined with 
CFHTLS-Deep photometric redshift catalogs for cluster 
detection. In closing, it must be emphasized that even 
though photometric redshift estimation may not be as 
expensive in integration time as spectroscopic redshifts, 
these methods still require deep, multi-band (> 5) imag- 
ing, while our method, K2, is designed to rely on just 
two colors and may thus find wider application. 

3. K2, A CLUSTER DETECTOR FOR MULTI-COLOR 
CFHTLS-W IMAGING 

Our detection method, K2, is designed to fully utilize 
the galaxy positions as well as their available colors in 
processed CFHTLS-W photometric catalogs. The prin- 
cipal idea is to identify simultaneous density enhance- 
ments in the projected 2D positions of the galaxies as 
well as in two colors independently; this method may 
be trivially extended to include more colors, if available. 
Cluster candidates are identified using a well defined met- 
ric computed with the position-color enhancement. Our 
search algorithm draws upon the strengths of four other 
methods reviewed in Section [21 with impro vements tai- 
lored for the CFHTLS-W data: -1- the C4 ([Miller et all 
120051) though without redshift informat ion, -2- to the 
Cut and Enhance (CE) (|Goto et al.ll200 2:) , with a galaxy 
based density enhan cement estimate, -3- the MaxBCG 
([Hansen et alJ I2005T ) without relying on a BCG evolu- 
tion model and -4- to the Cluster Red Sequence (CRS) 
([Gladders fc Yed 120001 ) though we stipulate a minimum 
of two colors to increase the robustness of the detections. 
Our method does not impose a luminosity distribution or 
a modeled radial distribution of cluster galaxies, thereby 
minimizing any selection bias in the detections. In ad- 
dition, we found that by locating the detection filter at 
each galaxy and not on a uniform grid defined over the 
survey area (as is done in both CE and in CRS meth- 
ods), our approach becomes more efficient in rejecting 
false positives due to noise peaks. This galaxy centered 
approach instead of a positional gri d has been adopted 
in both the Adaptiv e Matched Filter (i Kcpncr et al.lll99^ ) 
and the MaxBCG ([Hansen et al.l 120051 ) methods specifi- 
cally for this advantage. 

At the time of development of our cluster detection 
method, only gri photometry was available for most of 
CFHTLS-Wide survey area. To test the feasibility of ob- 
taining photometric redshifts with just these two colors, 
we compared SDSS spectroscopic redshifts (available for 
bright, r < 17 mag, galaxies in a fraction of the Wide 
fields) with photometr ic redshift estimates from HyperZ 
(Bolz onella et a l. 2000) using only the available gri mag- 
nitudes and errors. However, with only two colors the 
photometric redshift estimates for ~ 30% of the galaxies 
failed catastrophically. Hence K2 has been developed to 
work with just two color imaging and without reliance on 
redshift information; however, details of plans to incorpo- 
rate redshift information in our algorithm are discussed 
in Section [71 
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K2 is designed to work with processed photometric 
catalogs from multi-filter, wide field imaging. At present, 
we are using photometric catalogs for CFHTLS-W imag- 
ing which are available as part of the Terapix T05 data 
release (August 2008); Table [T] lists the numbers and 
other details of the available catalogs for the four survey 
fields of CFHTLS-W. The Terapix catalogs have b een 
produced with Sextractor (|Bertin fc Arnouti 119961 ) on 
the individual stacked image for each field in each fil- 
ter. During this processing, the i-image is used as the 
reference frame for source extraction from images in the 
other filters; the Terapix website 3 provides detailed ex- 
planations of all the image processing, astrometric and 
photometric calibrations during the production of the 
CFHTLS T05 catalogs. Since our cluster detection algo- 
rithm depends on the colors of galaxies, we independently 
compared the photometry of objects found in common 
between the CFHTLS-W fields and the SDSS imaging 
catalog; the differences between the photometry in these 
two catalogs were consistently smaller than the photo- 
metric errors, thus confirming that the Terapix photo- 
metric catalogs were suitable for our application. 

For each CFHTLS-W field, these Terapix gri photo- 
metric catalogs form the input to K2. For each object 
in the catalogs, we extract the ID number, sky positions 
(a, 6) in J2000 coordinates, x,y positions on the (20k x 
20k) image, the i-band half-light radius (HLR), and gri 
Sextractor Kron magnitudes with corresponding photo- 
metric errors; the Kron photometric apertures in all fil- 
ters are matched with those of the i-band image, which 
is used for object detection and segmentation. In order 
to exclude objects with poor photometry in one or more 
filters, we exclude all objects fainter than a user defined 
threshold or which have photometric errors greater then 
a user defined maximum; for CFHTLS-W, these selec- 
tion limits are i < 24 mag and magnitude errors < 0.5 
in all three filters, both limits being set to match the 
depth of the survey photometry. We then use the [i vs 
HLR) relation, a sample of which is shown in Figure [T] 
to isolate galaxies from stellar objects; the blue points, 
which lie below a HLR threshold demarcate the locus of 
point-like stellar objects. This HLR threshold, which is 
dependent on the PSF, is interactively selected for each 
image frame to account for field-to-field PSF variations; 
for the CFHTLS-W imaging, the HLR threshold is nom- 
inally ~ 0.6", while the upper limit is set at 4". 5. The 
magnitude limit for this star - galaxy discrimination us- 
ing the HLR method is also set by the user and is nom- 
inally i — 21.5mag; objects fainter than this limit are 
taken to be part of the galaxy population. Within the 
galaxy population, we select only a magnitude limited 
set of galaxies, referred as Bright Galaxies (BG), which 
we test for cluster membership; we use this magnitude 
selection to avoid saturated objects and those for which 
the photometric color errors are greater than the color 
cuts used for cluster member selection, as explained later 
in this section. For the CFHTLS-W, the bright end limit 
for the BG sample is set at i > 16 mag, while the faint 
end is i < 20 mag; as with the HLR limits, these magni- 
tude limits too may be varied interactively to take into 
account variations in the photometry between fields. The 
pertinent photometric and astrometric details of the BG 
are then written to a BG catalog for the cluster detection 
stage to follow. Similarly, a field catalog is prepared con- 



sisting of all objects in the observed field, which meet 
the magnitude and photometric selection criteria men- 
tioned above. Each 1 deg 2 CFHTLS- Wide field typically 
contains of order 2000 BGs and 150,000 field objects. 

For each BG in the catalog, K2 computes the (posi- 
tion+color) enhancement (referred as the weight) within 
a predefined aperture centered on the galaxy; the weight 
in each color is computed independently and the results 
are combined only in the last step of the detection pro- 
cess, when candidate cluster members are linked to iden- 
tify clusters. For each BG, we select all objects (field 
objects as well as other BGs) within a fixed aperture 
of angular size of 82", which corresponds to a radius 
0.5 Mpc at redshift z = 0.5; this physical radius cor- 
responds to half the virial radius of a typical cluster 
of virial mass, (M~ 10 14 M ). The chosen aperture 
size therefore adequately covers the cluster core from 
which (E+S0) galaxies contribute the majority of the po- 
sition and color wei g ht 8 . Other detect ion methods, e.g. 
IHansen et all (|2005D . IKim et al.l (|2002ft . adopt a similar 
detection aperture size for this reason. 

For the set of objects within the aperture, we then 
apply a color cut centered on the color of the BG to 
locate likely co-cluster members which lie along the Red 
Sequence. The width of the color cut, (= ±0.15 mag), 
corresponds to the photometric color error of a galaxy at 
the faint magnitude end (i — 24 mag) of our field sample. 
The cluster weight, as a measure of the BG being a cluster 
member (our null hypothesis), is computed as, 
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^ (5 d + e d ) (5 2 c + e c ) 



(1) 



where the position and color separations, 5, are measured 
with respect to the BG; suffixes refer to the separation 
in distance (d, measured in pixels) and in color (c, in 
magnitude). The summation is carried out only for the 
objects which lie within the aperture as well as within the 
imposed color separation from the BG; these position and 
color cuts normally yield a few tens of objects per BG. 
Softening parameters, e are used to prevent numerical 
errors due to division by zero; the softening parameter for 
distance, Cd — 0.02 pixel (a tenth of the astrometric error) 
and for color, e c = 0.005 mag (a tenth of the minimum 
color error for a BG ). The softening parameters thus 
represent a small fraction of the smallest expected error 
values in position and color. 

Next, we compute the field weight, which is a measure 
of the BG being just a field galaxy (the alternate hy- 
pothesis) . For this we use the field catalog with shuffled 
x, y positions for all the objects; their colors remain un- 
changed. The BG is placed at 1000 random positions 
in this shuffled field and at each position the weight, 
given by Equation [1] is computed. The median and the 
inter-quartile distance (IQD) of this set of 1000 weights 
is taken to be the field weight, Wf, and a measure of the 
scatter, o~f. Finally, the significance, which is the metric 
that we use to measure the likelihood of the BG being a 

8 We also tested the feasibility of basing aperture size on the 
photometric redshift of the BG; however, due to errors in the pho- 
tometric redshifts evaluated using only two colors, this adaptive 
aperture has not been implemented in the current version of K2. 
However, we aim to incorporate this approach in a future version 
tuned for five filter imaging, as discussed in [[7] 
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cluster member, is calculated by, 



w c - w f 



(2) 



We repeat this computation for all the galaxies in the 
BG catalog. 

Detected cluster members are those with significance 
values greater than a threshold, (3 a in the current im- 
plementation), in both the (g-r) and the {r-i) colors. 
This threshold value was determined using the detec- 
tion efficiency from our Monte-Carlo simulation results, 
described in H4.'2{ 

In the preceding steps, each BG is tested individually 
for cluster membership; subsets of these BGs may be 
members of the same clusters. We therefore group to- 
gether individual cluster members (BGs with significance 
greater t han the threshold) using a 'Friends-of-friends' al- 
gorithm (Huc hra fc Gell er 1982) in which galaxies within 
a certain linking length are combined into a single group. 
We use both colors and physical separations for linking 
member galaxies; our chosen link lengths are the widths 
of the color cuts and the radius of the detection aper- 
ture respectively. For clusters with multiple members, we 
identify the Bright Central Galaxy (BCG) as the mem- 
ber with the brightest i magnitude; for candidates in 
which only a single BG lies above the magnitude cut, 
that galaxy is also taken to be the BCG and its proper- 
ties are catalogued as described in S}6l Finally, we sepa- 
rate the detected clusters further into Gold, Silver and 
Bronze categories depending on their maximum detec- 
tion significances being greater than 5er in both colors, 
greater than 5 a in just one color or less than 5 a in both, 
respectively; note that all these cluster members satisfy 
the minimum detection significance of 3 tr, as mentioned 
earlier. 

3.1. Salient features of K2 

Before discussing the operational characteristics of K2, 
we highlight salient features which have been specifically 
designed to make the method flexible, portable and min- 
imally dependent on a fiducial cluster model. 

K2 has the built in flexibility to accept photometric 
catalogs in as many imaging colors as available, the min- 
imum being two; the greater the number of colors used, 
the higher is the overall detection significance, since only 
the combined significances (a minimum of 3er) in all col- 
ors flags a detection, as explained above. In addition, the 
results reported in this paper are from the current imple- 
mentation on the CFHTLS-W gri photometric catalogs. 
However, K2 is designed to accept photometric catalogs 
from any similar imaging survey, provided the necessary 
details {e.g. , HLR, sky positions, magnitudes) are avail- 
able. 

For ease of portability, our cluster detection routines 
are a group of IDL 9 procedures packaged into a sin- 
gle IDL project file. The complete cluster detection se- 
quence and preparation of cluster catalogs for each 1 deg 2 
CFHTLS-W field takes ^15 minutes on a single proces- 
sor PC. The only user intervention needed during exe- 
cution is to verify the star-galaxy discrimination step, 
similar to the one shown in Figure [T] where the default 

9 Integrated Development Language 



HLR and magnitude limits may be overridden to take 
into account field-to-field variations in the PSF and pho- 
tometry. 

It is to ensure that the dependency on a fiducial cluster 
model is kept to a minimum that we utilize the two-step 
detection process, testing only individual galaxies in Step 
1 to identify likely cluster members using the combined 
significances in two colors. The operational definition 
of a cluster member is an enhanced density of galax- 
ies within the detection aperture with colors matched 
(within the width of the color cut) to the color of the 
galaxy being tested. Only in Step 2 are the identified 
cluster member candidates then linked together, again 
using matched colors and the detection aperture as link- 
ing length, to identify member galaxies which belong to 
the same cluster. 

Therefore, K2 is capable of identifying cluster mem- 
bers with enhancements in any color and does not de- 
pend on the pres ence of a 'Red Sequ ence', as does 
the CRS method (|Gladders fc Yedl2000h . The method 
is thus capable of detecting even spiral rich struc- 
tures such as 'Hercule s cluster', Abell2215 at z=0.036, 
(jStruble fc Roodlll999h and the 'Blue Sequence' rep orted 
in the Virgo Cluster by Boselli and Gavazzi (2006). The 
detection significances of such clusters will consequently 
be higher only in blue colors, implying that (it — g) and 
(g — r) colors are better suited for the detection of such 
'blue' sequences. 

The above discussion may be extended to clusters in 
which a 'Red Sequence' is present but wit h a slope that is 
differe nt from the typical values giv en by IKodama et al.l 
(1998); an example is reported by lAdami et alj (|2007n 
from deep (B = 25.25 mag) observations of one quad- 
rant of the Coma cluster. In K2, the width of the color 
cut is set equal to a pre-factor x the fiducial slope of 
the CMR x the magnitude range used for detection; 
the pre-factor (default =1.5) is used to compensate for 
the associated photometric errors. The default color cut 
values implemented at present take into account typi- 
cal values of t h e CM R slope based on the results from 
IKodama et al.l (|1998h . It must be emphasized that the 
detection method does not fit the slope of the Red Se- 
quence, but merely identifies galaxies within the color 
cut. In order to identify clusters with CM R slopes with 
are ou tliers, such as the one reported by lAdami et al.l 
(|2007t ). an increased pre-factor and greater width of the 
color cut may be used. However, increasing the width 
of the color will consequently increase the contamination 
from foreground and background galaxies in the major- 
ity of clusters with typical CMR slopes and thus lower 
their detection significances. Fine tuning of this param- 
eter must therefore be assessed against the target cluster 
population and the intended application of the output 
catalog. 

A related comment applies to the effect of uncertain- 
ties in the measured magnitudes and colors of galaxies 
on the K2 detection efficiency. Since K2 computes the 
significance (Eq [2|) as a relative measure using cluster 
and field weights (both of which are affected equally by 
color errors), it is relatively insensitive to the method 
by which errors are estimated. For a given photometric 
catalog, the width of the color cut is chosen to match 
the color error of a galaxy at the faint end of the magni- 
tude range, as mentioned in Section [3l the default color 
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cut of ±0.3 mag applies to the current implementation 
on Terapix catalogs. 

Finally, K2 places a detection aperture on each BG 
and computes the weights, given by Eq [U using only 
the proximity of galaxies within the aperture and not 
by matching any radial profile. Since identified BGs are 
then linked in Step 2 using matched colors and the de- 
tection aperture as linking length, the overall effective 
aperture is much greater for multi-BG clusters than the 
single detection aperture. Thus, this Friends-of-friends 
linking process is capable of identifying clusters whose 
members are not distributed according to a classic pro- 
file, merging clusters with an extended distribution of 
members and of recovering filamentary structures such 
as the one observed in the i nfalling region in Abell 85 at 
z = 0.055 (Bouc ct al. 2008). If the colors of the member 
galaxies in such interacting clusters are different (greater 
than the width of the color cut), each is identified as an 
individual structure. Since K2 does not rely on a classic 
cluster profile, it identifies substructures within a cluster 
as part of the host cluster if their member galaxies' col- 
ors match, else as a separate entity provided the BGs in 
the substructure reach the detection threshold. We ver- 
ify these design features with Monte Carlo tests as well 
as using comparisons with published cluster catalogs in 
the following sections. 

4. MONTE CARLO TESTS FOR COMPLETENESS AND 
CONTAMINATION 

In the design of our detection method, we have kept 
the dependency on a cluster model to a minimum, thus 
decreasing bias against clusters which do not fit any pref- 
erential model. In addition, this low dependency on a 
cluster model minimizes the need for fine tuning of pa- 
rameters, and thus we are able to run the detection on 
all CFHTLS-W fields with the minimum of fine tuning. 

Despite these safeguards, the method is still prone to 
selection bias since it assumes the universal presence of 
the Red Sequence in all clusters, including poor clusters 
with few members. We test this assumption and thus es- 
timate the selection bias using two approaches: first with 
Monte Carlo tests using synthetic photometric catalogs 
as described in this section and then with a comparison 
with three other published CFHTLS cluster catalogs, the 
results of which are discussed in Section [5] In addition, 
we use these Monte Carlo simulations to estimate the 
contamination of our cluster catalog from false positives, 
e.g. chance superpositions of galaxies with matched col- 
ors, as well as completeness for a specified cluster richness 
and redshift. 

4.1. Simulation methodology 

For our Monte Carlo tests, we generated synthetic pho- 
tometric catalogs by placing synthetic clusters in a syn- 
thetic field of galaxies and stars, to mimic an observed 
CFHTLS-W field. With K2, the detection efficiency de- 
pends sensitively on color matches, hence we have en- 
sured that the synthetic BG and field photometric cata- 
logs are self-consistent in colors by using published and 
well tested empirical relations or theoretical results in 
constructing the synthetic catalogs. We also investi- 
gated the method using synthetic clusters on a shuffled 
CFHTLS-W field catalog and found that our recovery 
rate was unrealistically high (100% even at z ~ 1) be- 
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cause of significant differences between the synthetic col- 
ors of cluster members and the observed colors of the 
field galaxies. The following paragraph provides details 
of generating the synthetic field and BG catalogs used in 
our tests. 

For the synthetic field catalog, we assume a total field 
galaxy population of 150,000, the median number of 
galaxies with i < 24 mag found in a 1 deg 2 CFHTLS- 
W imaging field, see $3] We first generate the field 
galaxy number distribution, n(z, r'), as a function of red- 
shift and apparent r-magnitude, using the model source 
gal axy distribution comm only used in weak lensing stud- 
ies (|Brainerd et al.lll996f) 

This normalized, parametric fit to observed galaxy distri- 
butions in deep redshift surveys is applicable up to red- 
shift z ~ 1 and r-magnitude range, 20 < r < 24 mag, and 
is thus well matched to our range of inter est. We use the 
same functional form and values given in Brai nerd et al.l 
(1996) for the parameter, zq = zo(r) and the constant, 
P = 1.5. 

In this field galaxy population, we make the conser- 
vative assumption that there are equal numbers of early 
type E and SO and late type Sb galaxies. Observations 
of field galaxy luminosity functions and their evolution 
with redshift for different galaxy morphologies have how- 
ever shown that the field population is dominated by late 
types and that their luminosity functions evolve toward 
higher luminosity with redshift, while early and inter- 
medi ate types show little evolutionary trend (|Lin et al.l 
1999). Therefore our assumption of equal numbers as- 
signs a higher number density of (E+S0) galaxies to the 
field than the observed value; consequently, there is an 
increase in the field weight (more red galaxies in the field 
matching the color of the BGs) and therefore a reduc- 
tion in the detection significance; in other words, the de- 
tection significance from our Monte Carlo trials may be 
taken to be lower limits, making them a more stringent 
test of our detection method. 

For these three galaxy types, we compute the syn- 
thetic colors with the evolut ionary galaxy spectral ene rgy 
distribution (SED) models dBruzual fc Charlotll2003t) as- 
suming Scalo IMF (|Scaldll986l ). exponentially decaying 
star formation rates (with timescales of 1-2 Gyr for the 
early type galaxies and 5 Gyr for the Sb population) and 
a uniform formation redshift, z = 10. We have verified 
that varying the parameters of the SED models, such 
as formation redshift or metallicity, have little effect on 
the results of our simulations since such changes affect 
both the field and the cluster galaxy colors equally; for 
our simulations, it is necessary that the field and clus- 
ter galaxy colors be consistent relative to each other, 
their absolute values are of little import. These syn- 
thetic SEDs arc then convolved with the transmission 
functions of the CFHT gri filters. To the generated col- 
ors, we add photometric errors usin g the prescription in 
Sextractor (jBertin fc Arnoutsl 119961 ) with the published 
CFHTLS magnitude limits for the different filters. The 
positions of the field galaxies are distributed randomly 
within the (20kx20k) field. 

Our objective is to test the effectiveness of our method 
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in detecting cluster s with differ ent numbers of members 
and Abell richness (|Abel]Hl958ft and at various redshifts. 
Therefore, for generating each synthetic cluster, we start 
with a prescribed number of member galaxies and se- 
lect their absolute magnitudes from the cluster luminos- 
ity function. For our simulations, we use published val- 
ues of Schechter function parameters, M* — 5log(h) = 
—20.67, faint end slop e a = —1.2 and n ormalization, 
0* = 0.0146 h 3 Mpc~ 3 (jBlanton et al.l 12003ft . It must be 
mentioned that the Schechter parameters describing the 
cluster luminosity function have been observed to vary 
with the depth of the available imaging, e. g., a steeper 
faint end slope observed in deep surve ys dZucca et al. 
2006), as well as in detailed models (Khochfa r et al. 
2007f h To account for such variations, during the devel- 
opment of our Monte Carlo simulations we used several 
sets of Schechter parameters, (—20.2 < M* < —21.5, 
-0.7 < a < -1.25, 0.01 ^ 4>* < 0.045) taken from the 
literature (iPopesso et alj 120051 iZandivarez et al.l 120061 
iZucca et al.l 120061 iBlanton et al.1 120051 12003ft . However , 
our results for contamination and completeness of K2 
did not show any significant variations depending on the 
chosen values of the Schechter parameters. The results 
presented here therefore are obtained using only the sin- 
gle set of parameters given above. 

Since the objective of these Monte Carlo tests is to 
evaluate any inherent bias in our detection method only 
toward galaxy colors and not their absolute magnitudes, 
we take these field galaxy Schechter parameter values to 
apply equally well to our synthetic cluster member pop- 
ulation. With these generated absolute magnitudes, we 
then compute the Abell richness of the synthetic clus- 
ter at its assigned redshift. We iterate this procedure till 
the generated cluster richness matches the Abell richness 
value for which the simulation is carried out. 

We assign the morphologies and projected positions 
of member galaxies using observed results from liter- 
ature for the morphology-density relation and the lu- 
minosity segregation of cluster galaxies. Applying the 
observed luminosity segr egation from the ES Q Nearby 
Cluster Survey (ENACS. IBiviano et al.l (|2002ft ). we des- 
ignate all cluster members with absolute magnitude 
M r < —22 mag as elliptical galaxies, which popu- 
late only the cluster core. Probability density func- 
tions, derived from Figure 4 in IBiviano et al.l (|2002ft . are 
used to assign the projected radial distribution of these 
bright galaxies; their angular distribution is randomly 
assigned, assuming circular symmetry. I n keeping with 
these observational results provided by IBiviano et al.l 
(2002), cluster members fainter than M r = —22 mag 
are assumed to follow the Plummer radial distribution 
with a maximum radius, r max = 1.2 Mpc and core ra- 
dius, r c = 0.1r max ; their angular distribution is as- 
sumed to be circularly symmetric and is therefore ran- 
domly assigned. Depending on the radial position of 
each member, the morphology (E, SO and Sb) is de- 
tcrmi ned using the obser v ed mo rphology-density results 
from iThomas fc Katgertl (|2006f ) (the required probabil- 
ity densities are computed from results in Figure 2 of 
the publication). The colors of the E, SO and Sb clus- 
ter galaxies are obtained from the same SED models for 
these morphologies used for the field galaxies, thus main- 
taining color consistency. At the assumed cluster red- 
shift, member galaxies with apparent magnitudes below 



the limiting magnitude assumed in our work (i > 24 mag) 
are not included. Photometric errors for the cluster 
members are assigned using the same prescription as for 
the field galaxy population, thus maintaining consistency 
in their photometry. Finally, we wish to emphasize that 
though we use the Schechter luminosity function and 
the Plummer radial distribution function to generate the 
synthetic cluster, our detection method does not rely on 
these models and uses only cuts in position and in colors 
for locating candidate cluster members. 

The synthetic cluster is then merged with the field by 
placing it at a random position in the 20k x 20k field and 
converting the radial and angular positions of the cluster 
members into field x,y positions. All the cluster mem- 
bers brighter than the BG magnitude cut (see Section[3]) 
form the synthetic BG catalog, which is used for detec- 
tion; all the cluster members are treated as part of the 
field catalog as well. These two catalogs form the input 
to the cluster detector, thus following the methodology 
we use for the actual CFHTLS-W fields. The detection 
routine is repeated for a user defined number of trials 
(typically 100 times) for each synthetic cluster, with a 
different cluster location for each run. The (<? — r) and 
(r — i) detection significances from these trials as well 
as their median and inter quartile values are used for 
the analyses of the contamination rates and the com- 
pleteness of our detections These results from the Monte 
Carlo simulations are discussed in the following section. 

4.2. Monte Carlo results 

We first test the rate of false positives for different 
significance thresholds, (Equation^, in order to check 
for any underlying bias in our multi-color based cluster 
detection method, for instance due to a chance superpo- 
sition of galaxies being classified as a cluster. For this 
we select 100 random BGs with i < 20 mag and com- 
pute the median and scatter in just their field weights 
from 100 random locations in the field, repeating this 
process in redshift increments of Az = 0.05 in the in- 
terval 0.1 < z < 1. The false positive rate is taken 
to be the percentage number of times a galaxy in a ran- 
dom field position is flagged as a cluster with significance 
greater than the chosen threshold in a particular color. 
The histograms, shown separately for the (g — r) and the 
(r — i) colors in Figure [2) are compilations of the median 
false positive rates for a representative set of galaxies, 
each bar representing a redshift bin; results for different 
significance thresholds are represented by the colors in- 
dicated. We test the false detection rates for significance 
thresholds of 1, 2, 3 and 5 a. 

The results show that a 1 a detection threshold leads 
to an appreciable contamination due to false positives 
of ~ 20% in both colors; the contamination rate shows 
no clear correlation with redshift, since the detections are 
based on galaxy colors and not on magnitudes. The false 
positive rate decreases to ~ 5% at a detection threshold 
of 2 a and drops to 1 ~ 2% at 3 er; with the detection 
threshold set at 5 <r, the contamination drops essentially 
to zero. Based on these results, we have chosen a detec- 
tion threshold of 3 a. In addition, it must be emphasized 
that, for cluster detection, we impose the more stringent 
constraint that the significance should exceed the detec- 
tion threshold in both colors; therefore, the 3 a threshold 
in both colors reduces the combined contamination rate 
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to ~ 1%. It must be mentioned that in this False Posi- 
tives test, we included BGs with various magnitudes up 
to i < 20 mag and did not find any dependency. 

As a further test, we compute the completeness of our 
cluster detections using the recovery rate of synthetic 
clusters which are generated to resemble three observed 
clusters of different Abell richness classes. The results 
presented in Figure [3] pertain to clusters which resemble 
Coma (Abell richness class A c — 2), Fornax (A c = 1), 
and a poor WBL cluster (A c = 0) (similar to those in 
the catalog of poor clus ters by White, Blyton and Led- 
low. iWhite et al. (1999)); in the following discussion we 
refer to these three classes by the names of their tem- 
plate clusters. In the Plummer profile used for the radial 
distribution of the cluster members (see ^4. 1[) . the max- 
imum radius for each cluster class is scaled according to 
the richness, with Coma being assigned r max = 1.2 Mpc. 
The results presented here pertain to 100 realizations of 
each cluster class (ie. 100 Coma-like clusters generated 
using the method described in £14.11 and similarly for the 
other classes). The magnitudes and colors of each cluster 
as a function of redshift are computed taking into account 
cosmological dimming and k-corrections, as mentioned 
earlier; at each redshift, member galaxies fainter than 
the magnitude limit [i > 24 mag) , are not included in 
the computation. At each redshift, each cluster is placed 
at 100 random locations in the field and the detection 
significance is computed; the Monte Carlo significance 
for that cluster is taken to be the median of the values 
obtained from the 100 locations. This computation is 
repeated for each of the 100 realizations in each cluster 
class in redshift increments of Az — 0.05 in the interval 
0.1 < z < 1. 

Figure |4] plots the median significance values and the 
1 a scatter from the 100 realizations for each cluster class 
as a function of redshift (each point in that plot there- 
fore represents the median of 10 4 trials); the left and 
right panels represent the (g — r) and (r — i) colors used 
for detection. Over plotted in each panel is the 3 a de- 
tection threshold we have set using results from the False 
Positives tests. A comparison of the significances for the 
three classes at any given redshift shows the expected 
dominance of the Coma-like clusters with their larger 
population of (E+S0) galaxies, followed by the Fornax- 
like and WBL classes. As a function of redshift, the 
significance of each class trends downward as a greater 
number of cluster members drop below the magnitude 
threshold and therefore do not contribute to the com- 
puted significance. It is to be noted that the downward 
trend is not monotonic but shows spikes because the late 
type galaxies' colors fall within the BG color cut at these 
small redshift ranges, thus boosting the significance over 
that at even a slightly lower redshift. 

The completeness of each cluster class in each color 
at each redshift step is then computed as the fraction 
of realizations with median significance greater than 3, 
the threshold we have set for cluster detection based on 
the false detection estimates. The left and middle panels 
in Figure [3] show the completeness of our detections in 
the two colors independently. Since we stipulate that the 
significance must be greater than 3 a in both colors, the 
plots in rightmost panel show the combined completeness 
of our method for each of the three cluster classes as 



functions of redshift. 

The comparison shows that the detections are com- 
plete to > 80% for Coma in both the (g — r) and (r — i) 
colors at all redshifts up to z ~ 0.8; the detection effi- 
ciency drops rapidly at higher redshifts and asymptotes 
to a few percent (essentially equal to the false detection 
rate) by z ~ 0.9. This high detection efficiency is driven 
mainly by the bright early type galaxies in the core; their 
color clustering increases the detection significance while 
their bright magnitudes ensure that, even at the higher 
redshifts, they are above the i < 24 magnitude limit 
we have set for our detections. In addition, at all red- 
shifts, the (r — i) detection significance is higher, as is 
expected from the homogeneous red colors of the early 
type galaxies in the cluster core; the combined detection 
significance, (rightmost panel), therefore closely mimics 
the (5 — r) significance, which is the lower, therefore the 
more pertinent, of the two significances in determining 
the completeness. In the other cluster classes, the Fornax 
clusters are detected at efficiency > 80% only to z < 0.6 
while the WBL retain this efficiency only till z ~ 0.3. As 
pointed out earlier, the contribution of late type galaxies 
to the detection significance leads to the rise and fall in 
its trend seen in Figure IH instead of a monotonic de- 
crease. 

Our principal goal is to identify simultaneous enhance- 
ments of galaxy density and color as seen in photomet- 
ric data and use these as proxies to signal the pres- 
ence of mass enhancements due to clusters. With the 
Monte Carlo simulations, we have tested the efficacy of 
our multi-color based detection scheme and have also de- 
termined the optimal detection parameters for this ap- 
proach. 

5. COMPARISON WITH PUBLISHED CFHTLS CLUSTER 
CATALOGS 

As a further test of K2 as well as to understand any 
inherent detection bias, we compared our cluster cat- 
alogs for the CFHTLS-Dccp fields, which cover four 
patches each of 1 deg 2 within the larger Wide fields, with 
those of three other published cluster catalogs for the 
same sur vey regions, (1) the Matched Filter (MF) cluster 
catalogs (|Qlsen et all 120071 ) for the four CFHTLS-Deep 
fields, Dl - D4, (2) the X MM-LSS cluster ca talog for the 
ldeg 2 CFHTLS-D1 field (|Pierre et al.ll2007l) . and (3) the 
high likelihood cluster candidates i n the Adaptive Ker- 
nel (AK) catalogs for CFHTLS-D1 (jMazure et al.ll2007l ) 
based on th e photometric red shift catalogs for the four 
Deep fields (jllbert et all 120061 ). We discuss the compar- 
ison with each catalog in the f ollowing subsection s; it 
should be me n tioned that both lOlsen etlrtl (|2007l ) and 
iMazure et al.l (|2007j ) have carried out similar compar- 
isons of their optically selected catalogs with the X-ray 
selected XMM-LSS detections, thus providing a baseline 
against which we may compare our K2 candidates. This 
also helps us quantify the selection bias and completeness 
of these optical detection methods. 

5.1. Comparison with MF catalogs 

lOlsen et"al1 (2007) have produced their cluster catalogs 
by implementing t h e MF method, as originally proposed 
by iPostman et al.1 (|1996l ). to the 4 deg 2 covered by the 
four patches of CF HTLS-Deep field s. The filter uses the 
Schechter function (Schcchterl ll976f ) and a radial profile 
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chara cterized by a core and cutoff radii (Postman et al. 
1996) respectively for the luminosity and the radial dis- 
tribution of the cluster members in the fiducial cluster 
model; the values o f these model parameters are listed in 
IQlsen et al.1 ((2007) . They also describe the calibrations 
of the method using simulated data and provide details of 
their cluster catalogs for t he four CFHTLS-D eep fields. 
It must be mentioned that IQlsen et~a l. (2007) use visual 
inspection of RGB color images to classify the MF cluster 
candidates as follows: Class A candidates show a clear 
overdensity of galaxies with homogeneous colors; Class 
B show an overdensity but do not exhibit clustering in 
colors and, Class C visually indicate little enhancement 
eithe r in position or co lor of their likely member galax- 
ies. IQlsen et al.l ([2007) do not incorporate any of the 
available color information in the CFHTLS-D imaging 
in their detection algorithm, but only base the candidate 
cluster classification on their RGB colors. 

We summarize the MF detection results in Table [21 
classified according to their three categories. The num- 
bers listed under these three classes for the MF clusters 
are taken from their public catalogs 10 . Details of the 
numbers of clusters detected in each of the CFHTLS- 
Deep fields by our method, K2, are also given in Table 
[5] for comparison. The K2 detections are separated into 
Gold, Silver and Bronze categories (as defined in Section 

El). 

The results in Table [2] show that there are noticeable 
field to field variations both in the numbers of clusters 
detected by each method, as well as in the make-up of 
the categories of these detections, though the patterns of 
variations are different. Other than cosmic variance due 
the presence of large scale structure along the lines of 
sight of these 1 deg 2 patches, these differences may also 
arise due to variations in the depths of the observations 
in the three filters in these fields. 

Carrying this comparison a step further, Table [3] lists 
the breakdown of each detection class in the MF cluster 
catalog, which are also detected by our cluster method 
in this blind trial; the results pertain to the Deep 1 field, 
with similar breakdown values obtained for the other 
Deep fields. To generate this comparison, the positions of 
the candidate clusters in our cluster catalog are matched 
with those of the MF clusters, with a match being as- 
signed if the positions fall within a distance equal to our 
detection aperture. It must be noted that during this 
match up, it was found that ~ 30% of our clusters con- 
tain more than one MF cluster; this was specially true 
of our high significance Gold category. The reason is 
because we link together detected BGs within one de- 
tection aperture radius of each other into a single cluster 
candidate; the MF method on the other hand, treats each 
detection location as an independent cluster. 

The comparison shows that our detection method 
matches all the MF Class A detections that lie at red- 
shifts z < 0.8; the 7 non-detections lie at higher redshifts, 
according to the photometric rcdshift value assigned by 
the MF cluster filter. Similarly, the z > 0.8 MF candi- 
dates in Class B and C are not detected by our method 
along with ~ 25% of cluster detections at redshift be- 
low 0.8; nearly all these low redshift non-detections lie 
at z < 0.2. One possible cause may be that our detec- 

10 http://webviz.u-strasbg.fr/viz-bin/VizieR-2, J/A+A/461/81 



tion aperture size is not including all the galaxies in the 
core of these low redshift clusters to provide the 3 a min- 
imum detection significance. However, this effect is not 
observed in the Monte Carlo tests, where we observed 
the converse in that the detection significance was insen- 
sitive to aperture size above a certain minimum radius; 
this is the justification for the fixed detection aperture we 
use. An adaptive aperture based on photometric redshift, 
planned for a future version of cluster detector, may ad- 
dress this cause, if present. The other possible reason for 
the mismatch between MF and our detections at z < 0.2 
may be that the MF method, relying solely on galaxy 
density enhancement in a single filter, is more prone to 
contamination at lower redshifts, where the number of 
galaxies above the magnitude threshold is higher, lead- 
ing to enhanced density. 

Visual support of the above arguments is illustrated 
in Figure in which the three RGB images in the left 
column show cluster candidates detected by both the MF 
and K2 methods (top), by the MF only (middle) and by 
K2 only (bottom panel). The top panel clearly shows 
that a 'typical' cluster to the eye is detected by both 
methods consistently. On the other hand, a cluster with 
multiple central galaxies or without a radial symmetry 
(as seen in the lower panel) is missed by the MF perhaps 
due to the method's reliance on a fiducial cluster model; 
our detection method, which relies solely on position and 
color matching of the cluster members, is able to identify 
such cluster candidates. Our method, however, does not 
detect high redshift candidates (middle panel), mainly 
because of the faint magnitudes and associated color er- 
rors, which are greater than the color cuts used for select- 
ing member galaxies. These high redshift non-detections 
match the completeness trends predicted by our Monte 
Carlo results in Section 14. 2\ which show a steep drop in 
detection efficiency even for a rich Coma-like cluster at 
z > 0.8. In these high redshift clusters, the brighter z- 
magnitude of the cluster members may provide better 
S/N values, and therefore, the present redshift detection 
boundary of K2 may be extended to z > 1. with the 
future inclusion of the (i-z) color for detection. 

Finally, it must be emphasized that the MF uses only 
the j-filter photometric catalog during the detection pro- 
cess and thus has no color information, even though 
data for the Deep fields are available in all five filters; 
the only use of color in the MF method is to visually 
classify their detections into the three classes, as men- 
tioned above. Differences in detection efficiency between 
the cluster de tection me t hods have been shown to be 
significant by iKim et a l. (2002) in their comparison of 
the Matched Filter method with Vornoi Tessellation, the 
principal reason being the assumption of a uniform back- 
ground by the MF scheme. In our scheme, the inclusion 
of color makes it less susceptible to contamination by the 
background and this may be the reason for the differences 
in detection numbers as well as in the assignment of the 
detection classes. 

5.2. Comparison with XMM-LSS cluster catalogs 

The availability of the spectroscopically confirmed X- 
ray clusters from the XMM-LSS s urvey region whic h 
overlaps the CFHTLS-Deep 1 field [[Pierre et al.l [2007T) . 
provides a m ore stringent test f or our cluster detection 
method; since IQlsen et al.l (|2007| ) have carried out a sim- 
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ilar comparison with their Matched Filter catalogs, we 
also gain some insight into the detection bias of these 
two optical methods from these comparisons. We match 
the cluster detections in our catalogs with the published 
XMM-LSS cluster positions using the same method as 
for the MF comparison; a detected cluster BG in our 
catalog (with i < 20 mag, detection significances in both 
colors > 3 a) , whose separation from the published cen- 
tral position of X-ray emission (= position of the X-ray 
cluster) is less than the aperture radius used by K2 is 
taken to be a mat ched detection. We take these XMM- 
LSS positions from Olsen ct al. (2007) (Table 5) in order 
to maintain consistency in our comparison with the MF 
method. In preparing the cluster catalog for CFHTLS- 
Dl, we use the same K2 detection parameters as those 
used normally for all Wide fields and each cluster, if de- 
tected, is classified as Gold, Silver and Bronze based on 
their combined detection significances. The detection re- 
sults provided in Tabled] list whether each XMM cluster 
is detected by our method ( K2) and by the Mat ched Fil- 
ter (MF) scheme (as given in I Olsen et al.l (|2007f )) as well 
as the detection significances in the (g-r) and (r-i) colors 
and the corresponding cluster classification. As a visual 
check of this comparison, the three RGB images in the 
right column of Figure [5] show, (top panel) a represen- 
tative XMM-LSS cluster also detected by our method, 
(middle) an X-ray cluster not detected by K2 and (bot- 
tom panel) one of our Gold clusters in Deep-1, which 
does not have any associated X-ray emission and hence 
is not part of the XMM-LSS catalog. 

The comparison shows that our method detects 12 of 
the 17 XMM clusters. It is significant to note that our 
method detects all the X-ray clusters up to z — 0.8; all 
5 clusters which are not detected are faint, high redshift 
clusters at redshifts higher than 0.8. As seen in the mid- 
dle panel of Figure El the faint members of these high 
redshift clusters are indistinguishable against the field 
galaxies in the foreground; the i-magnitude cut used by 
K2 fails even to locate a BG for XLSSC029 at z = 1.05. 
On the other hand, the high significance Gold cluster in 
the bottom panel and which visually resembles a 'typical' 
cluster with a BCG and member galaxies with homoge- 
neous colors, is not part of the XMM-LSS catalog, per- 
haps due to a lack of associated X-ray emission. These re- 
sults are fully consistent with our results from the Monte- 
Carlo simulations in which the detection significances in 
both colors and therefore the completeness are > 80% 
even for Fornax-like clusters up to redshift 0.6 and then 
drop steeply off beyond z = 0.8 even for a Coma-like 
cluster with Abell richness = 2. However, as emphasized 
under Section 14.21 our principal focus at present, using 
only the available (g-r) and (r-i) colors, is to identify 
clusters up to redshift ~ 0.6; we aim to extend this to 
higher redshifts once we incorporate the z-filter catalogs 
into K2, as discussed in SJ7] 

In comparison, the MF success rate for the XMM 
clusters is 10 out of 17. It is interesting to note 
that the MF method fails to detect three clusters at 
z = 0.29, 0.34, and 0.46, all of which are detected with 
high significance by K2, while it successfully detects two 
higher redshift clusters at z = 0.92, and 1.05, though 
with redshift errors of ~ 0.1; at the same time the MF 
fails on a z = 0.92 cluster which is detected as a 'Gold' 
detection by our color based scheme. 



5.3. Comparison with Adaptive Kernel (AK) cluster 
catalogs 

iMazure et al.l (|2007f ) have used the publicly available 
photo metric redshift cat alogs for the four CFHTLS Deep 
fields (jllbert et al.ll2006f) to develop their cluster detec- 
tion method based on the Adaptive Kern el technique 
(|Silvermanlll986l . iDressler fc Shectmanlll98"8l ) and bench- 
mark their detections for the 1 deg 2 CFHTLS-D1 field 
against those i n the XMM-LSS cl uster catalogs for the 
same area. The lllbert et"aTI ()2006f ) catalogs provide well 
characterized photometric redshifts for all galaxies up to 
i = 24 mag in the four CFHTLS-Deep fields, obtained us- 
ing a combination of CFHTLS ugriz imaging, along with 
BVRI photometry fro m the VIMOS VLT Deep Survey 
(iLe Fevre et al . 2005) and infra-red J and JT-band imag- 
ing (jlovino et alJ "l2*005). The error in the photometric 
redshifts, estimated from a comparison with VVDS spec- 
troscopic redshifts for objects in the regions of overlap 
with the CFHTLS, is ±0.03 to 0.07 (depending on the 
galaxy magnitude) up to z = 1.5 . 

In using t he redshift catalog fo r CFHTLS-D1, for clus- 
ter detection lMazure et al.1 ([2007) split the galaxy sample 
into redshift bins, Az = 0.1 and apply the Adaptive Ker- 
nel to the resulting redshift slices to build density maps. 
Sextractor is then used to detect enhancements in galaxy 
density in these maps, as well as to obtain the geometric 
properties (ellipticity, orientation, etc.) of the regions of 
enhancement. In the 22 redshift slices used, a total of 40 
peaks are detected, which include all 8 Class 1 (CI) can - 
didates from the XMM-LSS catalogs ([Pierre et al.ll2006f ). 
Of these however, only 16 main structures appear in two 
or more redshift slices (including 4 XMM-LSS clusters) 
and are therefore identified as high likelihood cluster de- 
tections. Of these cluster candidates, 15 also have spec- 
troscopic redshifts by matching member galaxies with 
overlapping VVDS survey areas. 

Having already carried out an independent comparison 
of K2 detections with the XMM-LSS X-ray clusters, as 
described in £|5.2i we now compare our results with these 
16 high likelihood Adaptive Kernel detections; this pro- 
vides an additional benchmark for K2 against a different 
optical detection method based on photometric redshifts. 
It must be emphasized that the AK cluster catalog is built 
using high quality photometric redshifts, utilizing 8 col- 
ors plus sky positions for each object, against just the 2 
colors and positions on which K2 relies. 

The principal details of the AK candid ates are listed 
in Ta ble [5j reproduced from Table 6 in IMazure et al.l 
(2007); it must be pointed out that each Adaptive Kernel 
position given in the table refers to the position of the 
corresponding peak in the density map (determined by 
Sextractor) and not to an identified BCG of the candi- 
date cluster. Therefore, in order to compare detections 
with K2, which relies on galaxy positions, we compute 
the significance of each BG (refer to Eq [2] and Section 
[3] for details) within a co-moving aperture of IMpc at 
the redshift o f the c luster; this aperture size is used by 
IMazure et al.1 (|2007f ) in the comparison of their catalogs 
against the XMM-LSS, and we adopt the same aperture 
size for consistency, Cluster member candidates are then 
identified and linked by K2 using the methods explained 
in The last five columns in Table list whether K2 
detects the candidate, and if so, the significance values in 



12 



Thanjavur, Willis & Crampton 



gr and ri colors, the class of the candidate cluster and the 
positional offset (in arc sec) between the BCG identified 
by K2 and the position given in the AK catalog. 

The comparison shows that K2 detects 15 of the 16 
AK candidates up to the redshift, z ~ 0.85. This sig- 
nificant degree of match provides further support of the 
detection efficacy of K2 even given the limited color in- 
formation. At the same time however, this correspon- 
dence between the two detection methods is also to be 
expected given their reliance on the colors of the galax- 
ies, for detection in the case of K2 and for estimation of 
redshifts in the case of AK. In addition to the K2 detec- 
tions identified closest to the positions given in the AK 
catalogs, four apertures contain two K2 candidates each, 
while one aperture contains three K2 detections, These 
additional candidates have presumably not been segre- 
gated by AK as individual peaks due to the background 
threshold set in Sextractor or have not been linked by 
K2 due to scatter in their colors. The nature of these 
multiple detections may only be confirmed by spectro- 
scopic follow up. The only candidate not detected by K2 
is ID# 29 detected by AK in 2 redshift bins centered 
at z = 0.75 and 0.85; it must however be pointed out 
that the redshift of this candidate is unconfirmed since 
there is no overlap with the VVDS spectroscopic catalog. 
This uncertainty in the redshift of the candidate (likely 
a high redshift cluster) coupled with the rapid decline in 
the completeness of K2 with redshift (seen in Figure [3]), 
may be the reason for the non-detection by K2. These 
candidates may also be genuine sub-structures identified 
as a single cluster by the redshift based AK method, but 
segregated as individual structures due to differences in 
galaxy colors by K2. Finally, the AK candidate, ID# 12 
listed in Table which has no matching X-ray detec- 
tion in the XMM-LSS, is also detected by K2 as a Silver 
class candidate, as further support of the consistency of 
detections from these two different color based schemes. 

To summarize, the Monte Carlo results show that our 
cluster detection method is complete (> 80%) up to a 
redshift of 0.6 for clusters of richness matching or higher 
than Fornax; for poorer clusters, the completeness re- 
mains at > 80% until z=0.3 and declines to ~ 60% by 
z=0.6. These simulation results are supported by the 
comparison of our catalogs against clusters detected by 
the optical MF and AK methods as well as the XMM 
sample of X-ray selected clusters in the CFHTLS-D1. As 
discussed in SjTJ we expect to increase these completeness 
values even at higher redshifts with the inclusion of the 
(i — z) color in the detection process. 

6. CLUSTER CATALOGS FOR CFHTLS-W FIELDS 

We have constructed K2 cluster catalogs for all 
161 deg 2 of CFHTLS-W fields for which Terapix T05 
g, r and i photometric catalogs are available. Each 
1 deg 2 CFHTLS-W field has a corresponding K2 cata- 
log, named using the central RA and Dec of the field. 
In each catalog, we list the properties of the BCGs of 
the detected clusters and groups in that field; Section [3] 
explains how we identify the BCG of each object based 
upon the number of bright cluster members and their 
i-magnitudes. The properties listed in the K2 catalogs 
are the following: a unique ID number for each detection, 
J2000 sky positions (a and S), a photometric redshift es- 
timate (see Appendix for details), i magnitude and (g-r) 



and (r-i) colors of the BCG, the (g-r) and (r-i) detec- 
tion significances, the K2 detection class, the number of 
bright cluster members, the Abell richness and the num- 
ber of member galaxies, 72,32 as well as the limiting mag- 
nitude, TO32, used in computing the Abell richness. Of 
these, the positional and photometric details are taken 
directly from the Terapix catalogs, while all the other 
properties are estimated as part of the detection process. 
Table [6] shows the format of a sample K2 catalog. K2 
also lists the boundaries of each cluster candidate (in 
pixel coordinates on the 20k x 20k Megaprime image), 
which are used internally to generate color RGB images; 
we carry out the visual search for lenses on these RGB 
images 

In addition to these primary catalogs for all detections, 
K2 also generates an additional catalog for any cluster or 
group in which multiple BGs have been detected. These 
complementary catalogs list the astrometric and photo- 
metric details of the individual member galaxies as well 
as their estimated photometric redshifts. Table [7] shows 
a typical catalog of cluster member properties, of which 
the BCG properties form the first row. 

Figures [6] - [8] show comparisons of the numbers and 
characteristics of the detections in the four CFHTLS-W 
fields. Of these, Figure [5] shows the distribution of the 
K2 detection significances independently in the (g-r) and 
(r-i) colors. The histograms run from the minimum 3cr 
used as the threshold for detection to a maximum of 15 tr; 
this maximum value is chosen only for the purpose of the 
plot, and all objects with significances > 15 ct have been 
included in this bin. It must also be emphasized that an 
object is flagged as a detection only if it has the mini- 
mum 3 a significance in both colors; therefore, the (g-r) 
and (r-i) histograms for each CFHTLS-W region (corre- 
sponding left and right panels in Figure [5]) contain the 
same objects, though their significance distributions may 
differ. The annotated numbers on the right panels indi- 
cate the total number of detections in each region. Shown 
in Figure [7] are comparative histograms of the number of 
Gold, Silver and Bronze detections in each CFHTLS-W 
survey region; the histogram bins are in percentages of 
the total number of detections in that Wide field. Fig- 
ure [5] shows comparative histograms of the Abell richness 
values of the detections. The bins are labeled according 
to the generic names we have adopted for the three cate- 
gories, namely, group scale WBL (W), Fornax (F) and 
the rich Coma (C) class; the histograms are shown as 
percentages of the total number of detections in the cor- 
responding survey region. Finally, Table [51 provides a 
compilation of the quantitative and statistical compar- 
isons of the detections and their classifications in these 
four survey regions. In the following paragraphs, we dis- 
cuss these comparative characteristics of the candidates 
in the four survey regions. 

The histograms of detection significances (Figure \§§ 
show the expected trend of a larger number of candidates 
per significance bin at the lower end of the scale. How- 
ever, the histograms also illustrate that the total numbers 
of 5 a and greater detections are comparable to the num- 
bers of 3 — 5 a detections. This inference is reaffirmed 
by the comparison of detection classes in Figure [7J which 
clearly shows that there are equal numbers of high like- 
lihood Gold objects (^35% of the detections) as there 
are Bronze detections; as explained in £j3j the classifi- 
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cation scheme we use for these detections combines the 
significances in the (g-r) and (r-i) colors. 

We next compare the number of detections per square 
degree in each CFHTLS-W survey region. The median 
and inter-quartile distances, listed in Row 3 of Table [8j 
show that there are ~ 35 detections/ deg 2 (of all three 
classes) in the CFHTLS-W fields, and that the fields are 
statistically equivalent. The marginally lower value (30 
per deg 2 ) in W4 may be due to cosmic variance since 
the completeness limits in this region, given in Tabic [TJ 
match those of the others even though this region was 
added to the survey two years after its commencement., 
The small excess number of detections (~ 2 per deg 2 ) 
in the W3 region may also be due to the same reason, 
with a mass sheet lying along that line of sight. We note 
that 4 of the 9 lens candidates discovered so far in our 
search of gravitational lenses in the CFHTLS-W also lie 
in the W3 region, thus indicating an increased galaxy 
overdensity in that region. 

The median percentages in the three classes, shown in 
Figure [3 indicate that there are nearly equal numbers in 
the Gold and Bronze categories (^35% each), and are 
greater than the ~ 30% in the Silver class; this trend 
applies to all four Wide fields. This indicates that an 
object with a detection significance > 5 a in one color is 
likely also detected with a similar significance in the other 
color, and vice versa for the lower significance candidates. 

Finally, the distribution of the Abell richness classes 
of the candidates, Figure [H shows that the candidates 
are predominantly WBL-like objects. With the Abell 
richness as a proxy for mass, the statistics indicate that 
galaxy groups and poor clusters form the majority (> 
95%) of the total detections, or > 30 per deg 2 . This 
dominance is seen in all four Wide fields. Fornax and 
Coma-like rich clusters make up the remaining 5%, with 
a Fornax-class cluster found in every 1 — 2 deg 2 while 
a Coma-like cluster occurs every 2 — 4 deg 2 , with Wl 
showing the highest density (0.49 per deg 2 ). The total 
number of rich clusters in the K2 catalogs is 239, which 
provides us with a rich sample on which to focus our lens 
search. 

7. SUMMARY AND FUTURE DIRECTION 

We have developed a new cluster detection method, 
K2, specifically optimized for optical wide field, multi- 
color imaging. Based on the well tested Red Sequence 
approach, the algorithm uses the observed overdensity 
in optical colors and positions of early type galaxy popu- 
lation in the cores of galaxy groups and clusters to isolate 
them from the field galaxy population. The design draws 
upon select advantages of earlier cluster detection meth- 
ods with improvements tuned for the available imaging. 
By testing each galaxy for cluster membership, instead of 
running the detector on a uniform grid, false detections 
due to noise peaks or other imaging defects are mini- 
mized. Only galaxies in a magnitude limited sample are 
tested (the limit being set by the depth of photometry 
available) to improve execution speed and reduce false 
positives due to large photometric errors. 

K2 uses of a well defined metric to compute the com- 
bined position and color clustering; the metric, and the 
corresponding detection significance are calculated for 
each available color independently. Merging these inde- 
pendent lists in the final step of the process, and only 



selecting detections which lie above a well defined sig- 
nificance threshold, improves the likelihood of the can- 
didates being bona fide clusters. We define a consistent 
statistical estimate for the significance, and thus com- 
pute a quantitative measure of the likelihood of a galaxy 
being a cluster member. 

In order to provide the flexibility to readily incorpo- 
rate additional colors, K2 runs on the photometric cata- 
log for each color independently. Using additional colors, 
when available, improves the detection significance and 
the robustness of the detections at higher redshifts. The 
candidates are classified under Gold, Silver and Bronze 
categories to reflect their likelihood of being genuine clus- 
ters. The Abell richness and the photometric redshifts 
of the cluster candidates are computed as part of the 
detection process. 

Using Monte Carlo simulations, we have determined 
the optimal detection parameters for K2 as well as the 
completeness and contamination of the resulting clus- 
ter catalogs. Based on these tests, we determined that 
a fixed aperture for computing color and position en- 
hancements around each galaxy is more consistent than 
an adaptive one based on photometric redshift (this is 
because our tests indicated that when only two colors 
are used to estimate photometric redshifts of individual 
galaxies the failure rate is ~ 30%). The magnitude limit 
for the BG selection (i < 20mag), as well as the num- 
ber of trials (=1000) used to compute the field weight 
were also determined using these simulations. The Monte 
Carlo tests have shown that the contamination rate from 
false detections is less than 1% for the adopted minimum 
threshold of 3 o for the combined detection significances 
in both the (g — r) and (r — i) colors. The Monte-Carlo 
results have also shown that the completeness of our K2 
catalogs is > 80% for Fornax-like and richer clusters up 
to a redshift of 0.6, which is the redshift range of interest 
for searching for cluster lenses, our principal objective; 
for poor clusters of the WBL-class, the detection method 
is only complete at this level for redshifts up to 0.3 and 
drops to 60% by z=0.6. However, the lensing cross sec- 
tion of these low mass systems shows a corresponding 
steep decline with redshift and the likelihood of detect- 
ing lensed images in poor clusters at high redshift is low; 
therefore the lower completeness does not pose any dis- 
advantage for our application. 

As further verification, the K2 catalogs for the 
CFHTLS-Deep fields have been compared against three 
published cluster cata l ogs, n amely, the Matched Filter 
catalog by Ols en et ail (|2007[ ), the spectrosc opically con- 
firme d clusters in the XMM-LSS catalog (|Pierre et al.l 
2007), and cluster detections by t he Adaptive Kerne l 
method for the CFHTLS-D1 field (jMazure et al.ll200l 
which are based on photo metric redsh i ft cata logs for the 
CFHTLS-Deep fields by lllbert et al.l (poM ). In these 
comparisons, our candidates match all the confirmed X- 
ray clusters in the XMM-LSS up to a redshift of 0.8 as 
well as 15 of the 16 major structures detected by the 
Adaptive Kernel method in the same redshift range. At 
higher redshifts, the decrease in the number of matched 
detections is because the completeness of K2 drops off 
steeply above z ~ 0.8 due to two reasons: (1) the number 
of member galaxies which lie above the magnitude limit 
are too few to provide a detection significance above the 
established threshold value, and (2) in the current im- 
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plementation of K2 using only gri imaging, the Balmer 
break of the elliptical galaxies in the cluster cores, which 
are the principal contributors to the detection signifi- 
cance, falls outside wavelength range. To address this, 
we are in the process of incorporating additional (i — z) 
color in available fields and thus improving the complete- 
ness at redshifts above z ~ 0.8. Finally, the compari- 
son with the MF catalogs have also been consistent with 
the Monte Carlo estimates. Cases of mismatch arise due 
to the inherent selection biases present in the fiducial 
cluster-model based MF method and the color-based K2 
approach. 

Supported by these successful test results, we have 
generated K2 cluster catalogs for all the 161 deg 2 of 
CFHTLS-W fields for which Terapix T05 g, r and i-band 
photometric catalogs are available. Using an automated 
method, RGB color images of the high likelihood clus- 
ter candidates have been generated for a visual search 
for lensed arc images. We are also improving the photo- 
metric redshift values, at present determined using the 
2-color method described in the Appen dix, to an estimate 
using Hyper-Z (jBolzonella et al.ll2000r i with the full set 
of ugriz photometry in fields where it is available. The 
current version of K2 cluster catalogs are available on 
request from the authors. 
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APPENDIX 

PHOTOMETRIC REDSHIFTS OF K2 CLUSTER CANDIDATES 

By providing the photometric redshifts (photo-z) of the galaxy group and cluster candidates in our K2 catalogs, 
our objective is to increase their utility for other cosmological applications. Using the ava ilable two color photom etry, 
we initially attempted to determine the photo-z using available algorithms, e.g., HyperZ (Bol zonella et al.ll2000T l . and 
compared our estimates against the spectroscopic redshifts which were available in the SDSS spectroscopic catalogs 
for a small sample of our bright (r < 17.7mag) cluster member galaxies. Though the results were reasonable, our 
tests also indicated a significant failure rate of ~ 30%, mainly due to the uncertainty in determining the appropriate 
spectral type of the object using just two colors. 

Given this uncertainty due to limited color information, we have adopted a simpler, more direct approach to obtain 
an estimate of the photo-z of our candidates. Given that K2 detects clusters and groups by identifying member 
galaxies which lie along the Red Sequence, as described in Section [3J we make the reasonable assumption that all the 
member galaxies identified by K2 are early type E/SO galaxies. This assumption permits us to determine their photo-z 
directly using the evolution tracks of early type galaxy colors with redshift (refer to $2] for details and references). 
Specifically, we build and then u se a 2D fi tting function to estimate photo-z from the (g-r) and (r-i) colors of the 
member galaxies in our catalogs. Urn et all ()2002h have successfully used a similar approach for photo-z estimation in 
the Groth strip survey us ing just the one availa ble V — I color and I magnitude of the galaxies. The Cluster Red 
Sequence (CRS) method ( Glad ders k. Yedl2000f ) estimates redshifts based on this tight color-redshift correlation of 
cluster early type galaxies. In this Appendix, we describe our method, t he calibrations of our fitting function using 
the public release of photo-z catalogs for the four CFHTLS-Deep fields 11 (jllbert et al.ll2006l ) as well as error estimates 
using spectroscopic redshifts available in the public VIMOS VLT Deep Survey (VVDS) catalogs (|Le Fevre et al.ll2005h 
for the galaxy sample in the CFHTLS-D1 field (for brevity, we refer to the CFHTLS-Dccp photo-z catalogs as Ilbert 
catalogs and the spectroscopic ones as VVDS catalogs). The redshifts in the Ilbert catalogs are based on the five filter 
ugriz photometry available for the CFHTLS-Deep fields using their le Phare photo-z method (jllbert et al.ll2006h . 

We determined that a 2D cubic polynomial is adequate for tracing the evolution of the (g-r) and (r-i) colors of early 
type galaxies with redshift, as was shown by comparing the % 2 values of the fit with those for other fitting functions. 
In this fit, we restrict the redshift interval to z < 1, since the completeness of K2 drops off steeply above z = 0.8. 
The polynomial coefficients are obtained by fitting the available photo-z of early type galaxies in the Ilbert catalogs 
to their (g-r) and (r-i) colors, using a standard Levenberg-Marquardt least squares minimization technique. We are 
able to select this sample of early type galaxies since the galaxies in the Ilbert catalogs are flagged by spectral type, 
based upon the best fit template during photo-z estimation. In our fit for the coefficients, we only use galaxies where 
the confidence limit of the photometric redshift determination is listed as 95% or greater. For the CFHTLS-D1 field, 
this selection provides 3976 early type galaxies (with i < 24mag), with a similar number available for the other three 
fields. Each of the top two panels in Figure [9] plots (g-r) and (r-i) against the i-magnitude of the selected early type 
galaxies in the CFHTLS-D1; the different colors pertain to galaxies in equal Az — 0.2 intervals to highlight any redshift 

11 http://cencos.oamp.fr/ 
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dependent trends. The (r-i)-color shows a more systematic trend with i magnitude compared to the (g-r) color, which 
shows greater scatter mainly due to the degeneracy in this color at redshifts above 0.4 (as seen in the lower left panel). 
The lower two panels plot the photometric redshifts against each color, to highlight the redshift dependent trend to 
which we fit our chosen polynomial. The over-plotted points in these panels are the median and la scatter in the 
colors of galaxies binned by Az = 0.05 to better emphasize the trend with redshift. Also over-plotted as a solid line 
are the values of our fitting function evaluated at the median values in the bins (projected on to that color - redshift 
plane), to show the adequacy of the polynomial form for the fitting function. 

We use the spectroscopic redshifts in the public VVDS catalog for the CFHTLS-D1 field (jLe Fevrdl2007f l 12 to assess 
the errors in our photo-z determination. By matching (a, 5) positions, we identify our chosen sample of elliptical 
galaxies from the Ilbert catalog with the corresponding objects in the VVDS catalog; the VVDS catalog does not list 
their spectral type. Of the matched galaxies, we select only those which have been flagged with spectroscopic redshift 
confidence exceeding 95% and which lie below our chosen redshift limit of 1. The selection yields 363 galaxies in common 
between the two catalogs for the CFHTLS-D1 field. By comparing the spectroscopic and photometric redshifts of the 
selected early type galaxies, we determine that the median error in the photometric redshifts in the Ilbert catalog 
is 0.07, with no redshift dependency. It must be emphasized that the Ilbert photo-z values are determined using a 
combination of 8 optical and infra-red colors which are available for the CFHTLS-Deep fields (compared to the two 
colors in the CFHTLS-Wide fields, which we use for cluster detection and for the photo-z estimation). Therefore, the 
above error value provides a benchmark against which we next compare our redshift estimates as well as the associated 
errors for the same sample of galaxies . 

Using the (g-r) and (r-i) colors of the galaxies in the matched sample, we estimate their redshifts using our 2D 
polynomial. In this estimate, we use the combination of both colors because the slope of the (g-r) versus z correlation 
in the lower redshift range up to z ~ 0.4 is steeper than that of the (r-i) correlation (lower panels, Figure^ on the 
other hand, the (r-i) slope is steep at redshifts above 0.4 where the (g-r) color becomes degenerate. Therefore, using 
the colors in combination provides improved sensitivity up to z ~ 0.8, the redshift interval which our cluster candidates 
span. Figure [10] shows the correlation between our photo-z values and the corresponding spectroscopic redshifts from 
the VVDS catalog, (left panel) and the differences between the two estimates (right panel); over-plotted are also the 
median and standard deviations of the differences computed for galaxies in 0.2 redshift bins. For the sample of early 
type galaxies, the median error in our estimates is ~ 0.1, which is comparable to the median error in the Ilbert photo-z 
values for the same sample. However, given that these errors are estimated for a known set of early type galaxies, we 
take this error estimate of our method to be the lower limit since we do not take into account errors introduced by 
the uncertainty in the galaxy type and spectral template mismatch during photo-z determination. 

In applying our photo-z method to the cluster candidates in the K2 catalogs, we obtain the polynomial coefficients 
for each CFHTLS-Deep field independently to avoid any bias introduced due to color mismatch, e.g. from differing 
extinction values. We assume that these fits apply to the corresponding CFHTLS-W fields due to the overlap in sky 
coverage between the Deep and Wide surveys (i.e., the Dl fitting function is used for all the Wl fields, and similarly 
for the other fields). During the K2 cluster detection process, if a BG is flagged as a candidate cluster member, we 
obtain the photo-z for all the galaxies in the detection aperture which also match the BG colors (refer to §3] for details 
of the detection). The weighted mean and standard deviation of this set (with the «-band magnitudes as weights) are 
taken to be the photo-z of the BG and the associated error, thus improving the statistical confidence of our photo-z 
estimate (we also tested the alternative approach of sigma clipping to improve the photo-z estimate; however, as shown 
in Figure [S] a significant fraction of the detections are WBL-like poor clusters and groups with too few members to 
make this statistical approach feasible). The estimated photo-z value of the BG is written to the K2 catalogs along 
with the other properties. For single BG cluster candidates, this value is also taken to be the redshift of the cluster. For 
multi-BG candidates, we take the i-magnitude weighted mean and standard deviation of the photo-z of the member 
BGs to be the redshift of the host group or cluster. 
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Fig. 1. — Plot of i magnitude versus half light radius (HLR) of all objects in a 1 deg 2 CFHTLS-W field. The stellar locus defined by the 
HLR is used for star-galaxy segregation. The user defined magnitude and HLR selection cuts (described in Q are marked in dashed lines 
and the classification identifies stars (in blue), bright galaxies (red dots) and faint field galaxies (in green). Over- plotted are the magnitude 
cuts (chained lines) used for selecting the bright galaxies (BG), which are tested for cluster membership. 




Fig. 2. — Monte-Carlo estimates of the percentage of false positive detections as a function of redshift for different detection significance 
thresholds, as defined in Eq[2] The left and right panels refer to the (g — r) and (r — i) colors used for detection by K2. The plotted colors 
refer to threshold levels of la (red), 2<r (green), 3<r (blue) and 5<r (black, essentially zero) 
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Fig. 3. — Completeness as a function of redshift for a Coma-like (blue), Fornax-like (red) and a poor WBL cluster (black); the completeness 
values have been estimated for the (g-r) and (r-i) colors separately, as well as for both colors combined 
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Fig. 4. — The median and la scatter in the detection significance values as a function of redshift computed for the Coma-like, Fornax-like 
and a poor WBL cluster; the significance values are shown separately for the (g — r) and (r — i) colors. The red line indicates the 3cr 
detection threshold we have adopted for cluster membership 
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Fig. 5. — Mosaic of RGB color images of representative Matched Filter {left panels) and XMM-LSS clusters {right panels) matched 
against cluster candidates in our catalog for CFHTLS Deep-1 field. See jJ5]for a description of each panel. 




Fig. 6. — Histograms showing the distribution of the K2 detection significances of the cluster candidates in the (g — r) and (r — i) colors, 
shown independently, in the 161 deg 2 of CFHTLS-W imaging for which g,r and i imaging are available 




Fig. 7. — Histograms showing the relative percentages of Gold, Silver and Bronze cluster candidates in the four CFHTLS-W fields. The 
total number of candidates in each field is annotated for reference. 
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Fig. 8. — Histograms showing the relative percentages of Abell Richness values of the cluster candidates classified as WBL (W), Fornax 
(F) and Coma-like (C) clusters, corresponding to richness values -1, 0, 1 and higher, respectively. Annotated for comparison are the total 
number of candidates and areal coverage of each field along with the number of Fornax and Coma-like rich candidates, which have a higher 
likelihood of lensing background galaxies. 
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Fig. 9. — CMDs and correlations between photometric redshifts and colors of galaxies fitted with Elliptical galaxy templates for the 
CFHTLS-D1 field. The top panels show the CMDs for these galaxies, with colors representing redshift bins. The bottom panels are 
color versus photometric redshifts; the median and standard deviations of colors binned by 0.05 in redshift are over-plotted along with the 
projected fitted curve for each color (see Append ix lAl for details). All photometric and redshift values are taken from the public release of 
the photometric redshift catalogs by Ilbert et al. (2 0061 ). 
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Fig. 10. — Comparison of the redshifts obtained from our 3D fits with corresponding spectroscopic redshifts from the VVDS catalogs 
(left panel). The sample of 363 early type galaxies selected in the CFHTLS-D1 field is used for this comparison. The right panel shows 
the error in the fit ( = fit redshift - spectroscopic value) for the same set of galaxies; in order to check for any rcdshift dependency in this 
error, the median and lcr deviation values, estimated for galaxies in redshift bins of 0.2, are over-plotted. 



TABLE 1 

Completion statistics of CFHTLS-W survey 
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Completion statistics for the four CFHTLS-W fields computed using information provided at the Terapix website, http://terapix.iap.fr/. 
The columns listed are: N = number of fields observed; A = actual sky area covered, accounting for image boundaries; I = integration 
time; PSF = median seeing during observations; C s = Completeness limit for stellar objects; C g = Completeness limit for galaxies. The 



statistics quoted are the median values for each survey region in the corresponding filter, 
the K2 cluster catalogs are highlighted in bold font. 



The three filters used for the current version of 
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TABLE 2 

Comparison of cluster detections in CFHTLS-Deep fields by the Matched 
Filter, (MF) and K2 methods 
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Note. — Comparison of the numbers of clusters detected in the four CFHTLS - 
Deep fields by K2 with those in the Matched Filter Catalog (|Qlsen et al.H2007f ); 
for each detection method, the breakdown by detection classification as well as the 
total number are provided (see fJ3]and fJS]for details of the classification schemes 
used by each method) 



TABLE 3 

MF clusters in CFHTLS-D1 also detected by K2 



MF class 


Number 




-^"2 s H ver 


-^■^bronze 


K2 ND 


A class 


19 


12 






7 


B class 


13 


2 


3 


2 


6 


C class 


14 


6 


1 


1 


6 



Note. Breakdown of the MF cluster detections in 

CFHTLS-D1 field, which are also detected by our cluster 
detection method, K2, classified by the assigned class and as 
non-detections (ND); details of the classification scheme are 
described in Sj3] 



TABLE 4 

Comparison of K2 with XMM-LSS X-ray cluster detections in CFHTLS-D1 



XMM ID 


R. 


A. 


Dec. 




K2 


Sgr 




Class 


MF 


XLSSC029 


36 


.0172 


-4.2247 


1.05 


N 


0.00 


0.00 




N 


XLSSC044 


36 


.1410 


-4.2376 


0.26 


Y 


12.90 


9.38 


G 


Y 


XLSSJ022522. 7-042648 


36 


.3454 


-4.4468 


0.46 


Y 


20.52 


9.76 


G 


N 


XLSSC025 


36 


.3526 


-4.6791 


0.26 


Y 


9.27 


10.77 


G 


Y 


XLSSJ022529. 6-042547 


36 


.3733 


-4.4297 


0.92 


Y 


5.27 


9.23 


G 


N 


XLSSC041 


36 


.3777 


-4.2388 


0.14 


Y 


10.47 


10.58 


G 


Y 


XLSSC011 


36 


.5403 


-4.9684 


0.05 


Y 


11.08 


3.25 


S 


Y 


XLSSJ022609. 9-043120 


36 


.5421 


-4.5226 


0.82 


N 


2.51 


7.57 




N 


XLSSC017 


36 


.6174 


-4.9967 


0.38 


Y 


15.38 


10.13 


G 


Y 


XLSSC014 


36 


.6411 


-4.0633 


0.34 


Y 


14.27 


7.50 


G 


Y 


XLSSJ022651. 8-040956 


36 


.7164 


-4.1661 


0.34 


Y 


5.03 


3.37 


S 


N 


XLSSC005 


36 


.7877 


-4.3002 


1.05 


N 


-0.56 


-0.10 




Y 


XLSSC038 


36 


.8536 


-4.1920 


0.58 


Y 


3.26 


3.04 


B 


Y 


XLSSC013 


36 


.8588 


-4.5380 


0.31 


Y 


12.53 


9.66 


G 


Y 


XLSSC022 


36 


.9178 


-4.8586 


0.29 


Y 


14.39 


11.64 


G 


N 


XLSSJ022534.2-042535 


36 


.3925 


-4.4264 


0.92 


N 


1.01 


3.30 




Y 


XLSSC005b 


36 


.8000 


-4.2306 


1.00 


N 


1.07 


1.25 




N 



Note. - - Comparison of the detections of th e spectroscopically confirmed 
XMM-LSS X-ray clusters in the CFHTLS-D1 field (|Pierre et al.ll20OT ), using our 
detection method (K2) agains t the results with th e Matched Filter scheme (MF) 
for the same clusters given in Ols en et al.l ((2007). Our detections are classified 
further as Gold (G), Silver (S) or Bronze (B), as described in 33] 
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TABLE 5 

Comparison of Adaptive Kernel (Mazure et al. 200 7|) and K2 

CLUSTER DETECTIONS IN CFHTLS-D1 



AK ID 


R.A. 


Dec. 


z c 


K2 


Sgr 


Sri 


Class 


d 


1 


o u i 'von 

36.3789 


-4.2424 


0.138 


Y 


7.33 


5.97 


G 


91.84 


2 


Of? 7(101 

36. (981 


-4.1970 


0.185 


\r 

I 


9.37 


5.48 


c 

V T 


95.46 


3 


36.3/46 


-4.6831 


0.225 


V 
I 


5.00 


3.43 


t < 


75.55 


5 


36.6240 


-4.2523 


0.210 


\7 
1 


7.13 


5.11 


c 

V T 


81.21 


10 


36.3166 


-4.7515 


0.311 


Y 


13.15 


9.33 


G 


118.44 


6 


36.8416 


-4.5810 


0.308 


Y 


7.61 


4.55 


S 


51.65 


12 


36.6104 


-4.5286 


0.313 


Y 


5.32 


4.69 


S 


120.24 


15 


36.3842 


-4.2726 


0.542 


Y 


11.08 


5.88 


G 


37.46 


16 


36.8975 


-4.3768 


0.53 


Y 


5.71 


4.90 


S 


119.98 


17 


36.8481 


-4.6202 


0.543 


Y 


7.44 


5.36 


G 


101.74 


21 


36.4645 


-4.4997 


0.613 


Y 


5.97 


3.07 


S 


97.78 


19 


36.8646 


-4.5484 


0.610 


Y 


7.75 


5.84 


G 


66.33 


22 


36.6686 


-4.5096 


0.634 


Y 


4.28 


6.78 


S 


126.25 


28 


36.8791 


-4.2025 


0.784 


Y 


4.27 


7.11 


S 


94.77 


29 


36.5030 


-4.47142 




N 










34 


36.3925 


-4.4135 


0.920 


Y 


3.65 


4.42 


B 


90.72 



Note. — Comparison of the high likelihood cluster candidates in 
the AK catalogs with K2 detections in the CFHTLS-D1 field. The 
first four columns list the ID number, the central position in RA and 
Dec and the redshift of the AK de t ection s; these details are repro- 
duced from Table 6 in lMazure et al.| |2007j). The last five columns list 
whether K2 detects the candidate, and if so, the significance values in 
gr and ri colors, the class of the candidate cluster and the positional 
offset (in arc sec) between the BCG identified by K2 and the position 
given in the AK catalog. 



TABLE 6 

Format of master cluster catalog for each CFHTLS-W field 



ID 


a 


8 


2/« 


i 


g-r 


r — i 


Sgr 


S ,- i 


Class 


n cls 


7132 


A r 


m-32 


J140252+542504 


14:02:52.21 


+54:25:04.77 


0.279 


16.61 


0.887 


0.486 


4.64 


7.44 


1 


86 


83 


2 


19.27 


J140216+542445 


14:02:16.96 


+54:24:45.29 


0.221 


17.27 


0.736 


0.442 


4.95 


7.97 


1 


18 


16 


-1 


20.26 


J140356+542727 


14:03:56.55 


+54:27:27.26 


0.314 


16.69 


0.342 


0.238 


11.13 


12.22 


1 


36 


33 





19.35 


J140353+542728 


14:03:53.60 


+54:27:28.88 


0.257 


17.10 


0.892 


0.461 


9.00 


11.89 


1 


21 


29 


-1 


20.52 


J140646+542503 


14:06:46.01 


+54:25:03.16 


0.352 


17.51 


1.399 


0.614 


6.59 


5.10 


1 


3 


7 


-1 


21.27 



Note. - A sample listing of a typical K2 cluster catalog generated for each 1 deg 2 CFHTLS-W field. The 
properties listed are: a unique ID number, J2000 sky positions (a and S) , a photometric redshift estimate, i magnitude 
and (g-r) and (r-i) colors of the BCG, the (g-r) and (r-i) detection significances, the K2 detection class, the number of 
bright cluster members, the number of cluster members included in the richness class, the corresponding Abell richness 
and the limiting magnitude, 77132, used in computing the Abell richness. The positional and photometric details are 
taken directly from the Terapix catalogs; all the other properties are estimated during the detection process. Complete 
cluster catalogs for the 161 deg 2 CFHTLS-W fields are available on request from the authors 



K2: A new cluster detector 



TABLE 7 

Format of supplementary catalog for cluster member properties for each 

candidate cluster 



ID Class i g — rr — izfit x y HLR S gr S r i 

137343 1 16.607 0.887 0.486 0.271 16837.26 7955.87 10.64 12.64 7.44 

133713 1 17.444 0.914 0.565 0.294 16893.84 7718.00 7.53 4.42 4.65 

153526 1 19.021 0.719 0.568 0.269 16918.31 8337.70 7.09 9.22 8.45 

153427 1 19.038 0.897 0.565 0.287 16956.09 8362.39 5.90 7.25 7.59 



Note. — A sample listing of the supplementary catalog generated by K2 for any 
cluster or group in which multiple BGs have been detected. The columns listed are a 
unique ID number for each BG, K2 class of the cluster, i magnitude and (g-r) and (r-i) 
colors, photometric redshift, pixel positions in the corresponding Megaprime image 
and the {g-r) and (r-i) detection significances. Complete set of supplementary catalogs 
are available along with the corresponding master catalogs on request from the authors 



TABLE 8 

Detection statistics in K2 catalogs for each CFHTLS-W field 



Statistic 


Wl 


W2 


W3 


W4 


Number of fields 


72 


30 


49 


20 


Total detections 


2610 


1057 


1831 


646 


Detections/ dcg 2 


35 ± 12 


35 ± 9 


37 ±15 


30 ±14 


Class 


Gold 


984 


357 


713 


251 


Silver 


737 


301 


520 


163 


Bronze 


889 


399 


598 


232 


Gold [%] 


37 ± 16 


34 ± 13 


38 ±20 


40 ±16 


Silver [%] 


28 ± 16 


28 ± 16 


28 ± 19 


23 ±17 


Bronze [%] 


34 ± 18 


39 ± 16 


34 ±17 


37 ± 15 


Richness 


WBL total [%] 


2514 [96.3] 


1005 [95.1] 


1754 [95.8] 


632 [97.8] 


Fornax total [%] 


61 [2.3] 


44 [4.2] 


59 [3.2] 


14 [2.2] 


Coma total [%] 


35 [1.3] 


8[0.8] 


18 [1] 


0[-] 


WBL [number/ deg 2 ] 


34.9 


33.5 


35.796 


31.6 


Fornax [number/ deg 2 ] 


0.847 


1,467 


1.204 


0.7 


Coma [number/ dcg 2 ] 


0.486 


0.267 


0.367 




WBL [deg 2 /cluster] 


0.029 


0.03 


0.028 


0.032 


Fornax [ dcg 2 / cluster] 


1.18 


0.68 


0.83 


1.43 


Coma [dcg /cluster] 


2.06 


3.75 


2.72 





Note. Comparison of the detection statistics for the four 

CFHTLS-W fields deduced from the K2 cluster catalogs. The de- 
tails listed are: the total number of fields (= number of fields with 
Terapix T05 g, r and i photometric catalogs), total number of de- 
tections, number of detections/ deg 2 , detections in each detection 
class/ deg 2 , total number of detections in each Abell richness class, 
the surface density of each Abell richness class, and the inverse surface 
density (= sky area in which 1 cluster of that Abell richness class is 
found). All the quoted statistics are corresponding median and inter- 
quartile distances. The percentages listed in square brackets for the 
Abell richness classes (Rows 10 - 12) are computed using the total 
number of detections in that survey region, given on Row 2 of this 
table. 



